DirichletBC.cpp File Reference

#include <MoFEM.hpp>
#include <MethodForForceScaling.hpp>
#include <DirichletBC.hpp>

Go to the source code of this file.

Typedefs
using	FTensor1 = FTensor::Tensor1< double, 3 >

Functions
static MoFEMErrorCode	set_numered_dofs_on_ents (const Problem *problem_ptr, const FieldBitNumber bit_number, Range &ents, boost::function< MoFEMErrorCode(const boost::shared_ptr< MoFEM::NumeredDofEntity > &dof)> for_each_dof)
auto	get_rotation_from_vector (FTensor1 &t_omega)
auto	get_displacement (VectorDouble3 &coords, FTensor1 t_centr, FTensor1 t_normal, double theta)

Typedef Documentation

FTensor1

using FTensor1 = FTensor::Tensor1<double, 3>

Definition at line 90 of file DirichletBC.cpp.

Function Documentation

get_displacement()

auto get_displacement ( VectorDouble3 &  coords, FTensor1  t_centr, FTensor1  t_normal, double  theta  )

inline

Definition at line 115 of file DirichletBC.cpp.

                                                              {
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
 
  FTensor1 t_omega;
  FTensor1 t_coords(coords[0], coords[1], coords[2]);
  const double a = std::sqrt(t_normal(i) * t_normal(i));
  t_omega(i) = t_normal(i) * (theta / a);
  auto t_R = get_rotation_from_vector(t_omega);
  FTensor1 t_delta, t_disp;
  t_delta(i) = t_centr(i) - t_coords(i);
  t_disp(i) = t_delta(i) - t_R(i, j) * t_delta(j);
 
  return VectorDouble({t_disp(0), t_disp(1), t_disp(2)});
};

get_rotation_from_vector()

auto get_rotation_from_vector ( FTensor1 &  t_omega )

inline

Definition at line 92 of file DirichletBC.cpp.

                                                        {
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Tensor2<double, 3, 3> t_R;
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
  t_R(i, j) = t_kd(i, j);
 
  const double angle = std::sqrt(t_omega(i) * t_omega(i));
  if (std::abs(angle) <  1e-18)  
    return t_R;
 
  FTensor::Tensor2<double, 3, 3> t_Omega;
  t_Omega(i, j) = FTensor::levi_civita<double>(i, j, k) * t_omega(k);
  const double a = sin(angle) / angle;
  const double ss_2 = sin(angle / 2.);
  const double b = 2. * ss_2 * ss_2 / (angle * angle);
  t_R(i, j) += a * t_Omega(i, j);
  t_R(i, j) += b * t_Omega(i, k) * t_Omega(k, j);
 
  return t_R;
};

set_numered_dofs_on_ents()

static MoFEMErrorCode set_numered_dofs_on_ents ( const Problem *  problem_ptr, const FieldBitNumber  bit_number, Range &  ents, boost::function< MoFEMErrorCode(const boost::shared_ptr< MoFEM::NumeredDofEntity > &dof)>  for_each_dof  )

static

Definition at line 11 of file DirichletBC.cpp.

                      {
  MoFEMFunctionBegin;
 
  auto &dofs_by_uid = problem_ptr->getNumeredRowDofsPtr()->get<Unique_mi_tag>();
 
  for (auto eit = ents.pair_begin(); eit != ents.pair_end(); ++eit) {
 
    auto lo_dit = dofs_by_uid.lower_bound(
        DofEntity::getLoFieldEntityUId(bit_number, eit->first));
    auto hi_dit = dofs_by_uid.upper_bound(
        DofEntity::getHiFieldEntityUId(bit_number, eit->second));
 
    for (; lo_dit != hi_dit; ++lo_dit) {
      auto &dof = *lo_dit;
      if (dof->getHasLocalIndex())
        CHKERR for_each_dof(dof);
    }
  }
 
  MoFEMFunctionReturn(0);
};