ElasticTie::OpTiePrismRhs Struct Reference

#include "tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp"

Inheritance diagram for ElasticTie::OpTiePrismRhs:

Collaboration diagram for ElasticTie::OpTiePrismRhs:

Public Member Functions
	OpTiePrismRhs (boost::shared_ptr< MatrixDouble > disp_ptr)
double	getSlaveFaceArea ()

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data) override

Private Attributes
boost::shared_ptr< MatrixDouble >	dispPtr

Detailed Description

Definition at line 563 of file ElasticTie.hpp.

Constructor & Destructor Documentation

OpTiePrismRhs()

ElasticTie::OpTiePrismRhs::OpTiePrismRhs ( boost::shared_ptr< MatrixDouble >  disp_ptr )

inline

Definition at line 565 of file ElasticTie.hpp.

      : TiePrismRhsBase("LAMBDA", "LAMBDA", OPROW, nullptr),
        dispPtr(disp_ptr) {
    if (!dispPtr)
      THROW_MESSAGE("Pointers are not set");
  }

Member Function Documentation

getSlaveFaceArea()

double ElasticTie::OpTiePrismRhs::getSlaveFaceArea ( )

inline

Definition at line 572 of file ElasticTie.hpp.

{ return getAreaF3(); }

iNtegrate()

MoFEMErrorCode ElasticTie::OpTiePrismRhs::iNtegrate ( EntitiesFieldData::EntData &  row_data )

inlineoverrideprotected

Definition at line 575 of file ElasticTie.hpp.

                                                                        {
    MoFEMFunctionBeginHot;
 
    const auto *fe =
        dynamic_cast<const FlatPrismElementForcesAndSourcesCore *>(ptrFE);
    if (!fe)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "TIE prism RHS operator is attached to a wrong finite element");
 
    const auto nb_gauss_pts = getGaussPts().size2();
    if (!nb_gauss_pts || nb_gauss_pts % 2) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Expected concatenated [slave|master] TIE Gauss points");
    }
    const auto nb_slave_gauss = nb_gauss_pts / 2;
    if (dispPtr->size1() != nb_gauss_pts) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "TIE displacement values at Gauss points were not prepared "
              "before RHS assembly");
    }
 
    if (nbRows % SPACE_DIM)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected number of LAMBDA dofs on TIE prism rows");
 
    auto &row_n = row_data.getN();
    if (!row_n.size1() || !row_n.size2())
      MoFEMFunctionReturnHot(0);
    if (row_n.size1() != nb_gauss_pts)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected number of Gauss points on TIE multiplier block");
 
    const int nb_row_base = nbRows / SPACE_DIM;
    if (static_cast<size_t>(nb_row_base) > row_n.size2())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "More TIE multiplier bases than shape functions");
    if (dispPtr->size2() != SPACE_DIM)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected TIE displacement matrix size: size2=%d expected=%d",
              static_cast<int>(dispPtr->size2()), SPACE_DIM);
    if (locF.size() != static_cast<size_t>(nbRows))
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected TIE RHS block size: rowSide=%d rowType=%d nbRows=%d "
              "locF.size()=%d row_n=(%d,%d) disp=(%d,%d)",
              rowSide, rowType, nbRows, static_cast<int>(locF.size()),
              static_cast<int>(row_n.size1()), static_cast<int>(row_n.size2()),
              static_cast<int>(dispPtr->size1()),
              static_cast<int>(dispPtr->size2()));
 
    auto t_w = getFTensor0IntegrationWeight();
    for (size_t gg = 0; gg != nb_slave_gauss; ++gg) {
      const double slave_alpha = getSlaveFaceArea() * t_w;
      auto t_row_base = row_data.getFTensor0N(gg, 0);
      for (int rr = 0; rr != nb_row_base; ++rr) {
        const double row_base = t_row_base;
        if (row_base != 0.) {
          auto t_row_rhs = getFTensor1FromPtr<SPACE_DIM>(&locF[SPACE_DIM * rr]);
          for (int dd = 0; dd != SPACE_DIM; ++dd) {
            t_row_rhs(dd) += row_base * slave_alpha *
                             ((*dispPtr)(gg, dd) -
                              (*dispPtr)(gg + nb_slave_gauss, dd));
          }
        }
        ++t_row_base;
      }
 
      ++t_w;
    }
 
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

dispPtr

boost::shared_ptr<MatrixDouble> ElasticTie::OpTiePrismRhs::dispPtr

private

Definition at line 648 of file ElasticTie.hpp.

---

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

• tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp