OpElasticTools::OpEssentialRHS_U Struct Reference

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

Inheritance diagram for OpElasticTools::OpEssentialRHS_U:

Collaboration diagram for OpElasticTools::OpEssentialRHS_U:

Public Member Functions
	OpEssentialRHS_U (std::string disp_name, boost::shared_ptr< CommonData > common_data_ptr, VectorDouble &disp, const Range &essential_ents)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
VectorDouble &	disP
Range	essentialEnts

Detailed Description

Definition at line 949 of file ElasticOperators.hpp.

Constructor & Destructor Documentation

OpEssentialRHS_U()

OpElasticTools::OpEssentialRHS_U::OpEssentialRHS_U ( std::string  disp_name, boost::shared_ptr< CommonData >  common_data_ptr, VectorDouble &  disp, const Range &  essential_ents  )

inline

Definition at line 951 of file ElasticOperators.hpp.

      : BoundaryEleOpAssembly(disp_name, disp_name, BoundaryEleOp::OPROW),
        essentialEnts(essential_ents), disP(disp),
        commonDataPtr(common_data_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpElasticTools::OpEssentialRHS_U::iNtegrate ( EntitiesFieldData::EntData &  row_data )

inline

Definition at line 958 of file ElasticOperators.hpp.

                                                               {
    MoFEMFunctionBegin;
    const int nb_row_dofs = row_data.getIndices().size();
    const int nb_row_base_functions = row_data.getN().size2();
 
    if (nb_row_dofs) {
 
      EntityHandle ent = getFEEntityHandle();
      if (essentialEnts.find(ent) == essentialEnts.end()) {
        MoFEMFunctionReturnHot(0);
      }
 
      auto t_disp = getFTensor1FromMat<3>(*commonDataPtr->mDispPtr);
 
      const int nb_integration_pts = getGaussPts().size2();
      auto t_row_val = row_data.getFTensor0N();
      auto t_w = getFTensor0IntegrationWeight();
      const double vol = getMeasure();
      auto t_coords = getFTensor1CoordsAtGaussPts();
 
      Tensor1<double, 3> bc_disp(disP(0), disP(1), disP(2));
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        const double a = vol * t_w;
        
        // // for testing radial displacement
        // auto theta = atan2(t_coords(1), t_coords(0));
        // bc_disp(0) = disP(0) * cos(theta);
        // bc_disp(1) = disP(1) * sin(theta);
 
        auto t_nf = getFTensor1FromArray<3, 3>(locF);
        int rr = 0; 
        for (;rr != nb_row_dofs / 3; ++rr) {
          t_nf(i) += a * t_row_val * (t_disp(i) - bc_disp(i));
          ++t_row_val;
          ++t_nf;
        }
        for (; rr < nb_row_base_functions; ++rr)
          ++t_row_val;
        ++t_disp;
        ++t_coords;
        ++t_w;
      }
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpElasticTools::OpEssentialRHS_U::commonDataPtr

private

Definition at line 1005 of file ElasticOperators.hpp.

disP

VectorDouble& OpElasticTools::OpEssentialRHS_U::disP

private

Definition at line 1006 of file ElasticOperators.hpp.

essentialEnts

Range OpElasticTools::OpEssentialRHS_U::essentialEnts

private

Definition at line 1007 of file ElasticOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/ElasticOperators.hpp