OpBrokenTractionBc Struct Reference

#include <users_modules/eshelbian_plasticit/src/EshelbianOperators.hpp>

Inheritance diagram for OpBrokenTractionBc:

Collaboration diagram for OpBrokenTractionBc:

Public Member Functions
	OpBrokenTractionBc (const std::string row_field, boost::shared_ptr< TractionBcVec > bc_data, boost::shared_ptr< double > piola_scale_ptr, std::vector< boost::shared_ptr< ScalingMethod >> smv)
MoFEMErrorCode	iNtegrate (EntData &data)

Protected Attributes
boost::shared_ptr< TractionBcVec >	bcData
boost::shared_ptr< double >	piolaScalePtr
std::vector< boost::shared_ptr< ScalingMethod > >	scalingMethodsVec

Detailed Description

Examples

EshelbianPlasticity.cpp.

Definition at line 334 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpBrokenTractionBc()

OpBrokenTractionBc::OpBrokenTractionBc ( const std::string  row_field, boost::shared_ptr< TractionBcVec >  bc_data, boost::shared_ptr< double >  piola_scale_ptr, std::vector< boost::shared_ptr< ScalingMethod >>  smv  )

inline

Definition at line 336 of file EshelbianOperators.hpp.

      : FormsIntegrators<FaceUserDataOperator>::Assembly<PETSC>::OpBase(
            row_field, row_field, FaceUserDataOperator::OPROW),
        bcData(bc_data), piolaScalePtr(), scalingMethodsVec(smv) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpBrokenTractionBc::iNtegrate

(

EntData &

data

)

Examples

EshelbianOperators.cpp.

Definition at line 1278 of file EshelbianOperators.cpp.

                                                          {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(3, i);
 
  int nb_dofs = data.getFieldData().size();
  int nb_integration_pts = getGaussPts().size2();
  int nb_base_functions = data.getN().size2();
 
  double ts_t = getFEMethod()->ts_t;
  if (EshelbianCore::dynamicRelaxation) {
    ts_t = EshelbianCore::dynamicTime;
  }
 
  double time_scale = 1;
  for (auto &sm : scalingMethodsVec) {
    time_scale *= sm->getScale(ts_t);
  }
 
#ifndef NDEBUG
  if (this->locF.size() != nb_dofs)
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Size of locF %d != nb_dofs %d", this->locF.size(), nb_dofs);
#endif // NDEBUG
 
  auto integrate_rhs = [&](auto &bc, auto calc_tau) {
    MoFEMFunctionBegin;
 
    auto t_val = getFTensor1FromPtr<3>(&*bc.vals.begin());
    auto t_row_base = data.getFTensor0N();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    double scale = (piolaScalePtr) ? 1. / (*piolaScalePtr) : 1.0;
 
    for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
      const auto tau = calc_tau(t_coords(0), t_coords(1), t_coords(2));
      auto t_f = getFTensor1FromPtr<3>(&*this->locF.begin());
      int rr = 0;
      for (; rr != nb_dofs / SPACE_DIM; ++rr) {
        t_f(i) -= time_scale * t_w * t_row_base * tau * (t_val(i) * scale);
        ++t_row_base;
        ++t_f;
      }
 
      for (; rr != nb_base_functions; ++rr)
        ++t_row_base;
      ++t_w;
      ++t_coords;
    }
    // Multiply by 2 since we integrate on triangle, and hybrid constrain is
    // multiplied by norm.
    this->locF *= 2. * getMeasure();
    MoFEMFunctionReturn(0);
  };
 
  // get entity of face
  EntityHandle fe_ent = getFEEntityHandle();
  for (auto &bc : *(bcData)) {
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
 
      int nb_dofs = data.getFieldData().size();
      if (nb_dofs) {
 
        if (std::regex_match(bc.blockName, std::regex(".*COOK.*"))) {
          auto calc_tau = [](double, double y, double) {
            y -= 44;
            y /= (60 - 44);
            return -y * (y - 1) / 0.25;
          };
          CHKERR integrate_rhs(bc, calc_tau);
        } else {
          CHKERR integrate_rhs(bc, [](double, double, double) { return 1; });
        }
      }
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcData

boost::shared_ptr<TractionBcVec> OpBrokenTractionBc::bcData

protected

Definition at line 347 of file EshelbianOperators.hpp.

piolaScalePtr

boost::shared_ptr<double> OpBrokenTractionBc::piolaScalePtr

protected

Definition at line 348 of file EshelbianOperators.hpp.

scalingMethodsVec

std::vector<boost::shared_ptr<ScalingMethod> > OpBrokenTractionBc::scalingMethodsVec

protected

Definition at line 349 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticit/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticit/src/impl/EshelbianOperators.cpp