OpRotatingFrameTools::OpRotatingFrameBoundaryRhs Struct Reference

#include <users_modules/multifield_plasticity/src/RotatingFrameOperators.hpp>

Inheritance diagram for OpRotatingFrameTools::OpRotatingFrameBoundaryRhs:

Collaboration diagram for OpRotatingFrameTools::OpRotatingFrameBoundaryRhs:

Public Member Functions
	OpRotatingFrameBoundaryRhs (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< DomainSideEle > side_op_fe)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
boost::shared_ptr< DomainSideEle >	sideOpFe

Detailed Description

Definition at line 29 of file RotatingFrameOperators.hpp.

Constructor & Destructor Documentation

OpRotatingFrameBoundaryRhs()

OpRotatingFrameTools::OpRotatingFrameBoundaryRhs::OpRotatingFrameBoundaryRhs	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< DomainSideEle >	side_op_fe
	)

Definition at line 162 of file RotatingFrameOperators.cpp.

    : BoundaryEleOpAssembly(field_name, field_name, OPROW),
      commonDataPtr(common_data_ptr), sideOpFe(side_op_fe) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpRotatingFrameTools::OpRotatingFrameBoundaryRhs::iNtegrate

(

EntitiesFieldData::EntData &

data

)

Definition at line 168 of file RotatingFrameOperators.cpp.

                                                                  {
  MoFEMFunctionBegin;
 
  const size_t nb_gauss_pts = getGaussPts().size2();
  const size_t nb_dofs = data.getIndices().size();
 
  if (nb_dofs) {
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    CHKERR loopSideVolumes("dFE", *sideOpFe);
    auto t_grad = getFTensor2FromMat<3, 3>(*(commonDataPtr->mGradPtr));
    size_t nb_base_functions = data.getN().size2();
    auto t_base = data.getFTensor0N();
    auto t_normal = getFTensor1Normal();
    const double ls = sqrt(t_normal(i) * t_normal(i));
    t_normal(i) /= ls;
 
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
      auto t_nf = getFTensor1FromArray<3, 3>(locF);
 
      if (getNormalsAtGaussPts().size1() == nb_gauss_pts) {
        VectorDouble n = getNormalsAtGaussPts(gg);
        auto t_n = getFTensor1FromPtr<3>(&*n.data().begin());
        t_normal(i) = t_n(i) / sqrt(t_n(j) * t_n(j));
      }
 
      const double alpha = t_w * getMeasure() * (*cache).density;
      Tensor1<double, 3> omg;
      omg(i) = (*cache).Omega(i, j) * t_coords(j);
 
      size_t bb = 0;
      for (; bb != nb_dofs / 3; ++bb) {
        // cerr << test << endl;
        t_nf(i) += alpha * t_base *
                   ((t_grad(i, j) + kronecker_delta(i, j)) * omg(j)) *
                   (omg(k) * t_normal(k));
        // alpha * t_base * (t_grad(i, j) * omg(j));
 
        ++t_nf;
        ++t_base;
      }
      for (; bb < nb_base_functions; ++bb)
        ++t_base;
 
      ++t_grad;
      ++t_coords;
      ++t_w;
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpRotatingFrameTools::OpRotatingFrameBoundaryRhs::commonDataPtr

private

Definition at line 37 of file RotatingFrameOperators.hpp.

sideOpFe

boost::shared_ptr<DomainSideEle> OpRotatingFrameTools::OpRotatingFrameBoundaryRhs::sideOpFe

private

Definition at line 38 of file RotatingFrameOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/RotatingFrameOperators.hpp
• users_modules/multifield_plasticity/src/impl/RotatingFrameOperators.cpp