EshelbianPlasticity::OpBrokenAnalyticalTractionBc Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianPlasticity.hpp"

Inheritance diagram for EshelbianPlasticity::OpBrokenAnalyticalTractionBc:

Collaboration diagram for EshelbianPlasticity::OpBrokenAnalyticalTractionBc:

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

Protected Attributes
boost::shared_ptr< AnalyticalTractionBcVec >	bcData
boost::shared_ptr< double >	piolaScalePtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap

Detailed Description

Definition at line 559 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

OpBrokenAnalyticalTractionBc()

EshelbianPlasticity::OpBrokenAnalyticalTractionBc::OpBrokenAnalyticalTractionBc ( const std::string row_field, boost::shared_ptr< AnalyticalTractionBcVec > bc_data, boost::shared_ptr< double > piola_scale_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv )

inline

Definition at line 561 of file EshelbianPlasticity.hpp.

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpBrokenAnalyticalTractionBc::iNtegrate

(

EntData &

data

)

Examples

EshelbianOperators.cpp.

Definition at line 2205 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 time = getFEMethod()->ts_t;
  if (EshelbianCore::dynamicRelaxation) {
    time = EshelbianCore::dynamicTime;
  }
 
#ifndef NDEBUG
  if (this->locF.size() != nb_dofs)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Size of locF %ld != nb_dofs %d", this->locF.size(), nb_dofs);
#endif // NDEBUG
 
  // get entity of face
  EntityHandle fe_ent = getFEEntityHandle();
  for (auto &bc : *(bcData)) {
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
 
      MatrixDouble analytical_expr;
      // placeholder to pass boundary block id to python
      auto [block_name, v_analytical_expr] =
          getAnalyticalExpr(this, analytical_expr, bc.blockName);
      auto t_val =
          getFTensor1FromMat<3>(v_analytical_expr);
      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) {
 
        auto t_f = getFTensor1FromPtr<3>(&*this->locF.begin());
        int rr = 0;
        for (; rr != nb_dofs / SPACE_DIM; ++rr) {
          t_f(i) -= t_w * t_row_base * (t_val(i) * scale);
          ++t_row_base;
          ++t_f;
        }
 
        for (; rr != nb_base_functions; ++rr)
          ++t_row_base;
        ++t_w;
        ++t_coords;
        ++t_val;
      }
      this->locF *= getMeasure();
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcData

boost::shared_ptr<AnalyticalTractionBcVec> EshelbianPlasticity::OpBrokenAnalyticalTractionBc::bcData

protected

Examples

EshelbianOperators.cpp.

Definition at line 572 of file EshelbianPlasticity.hpp.

piolaScalePtr

boost::shared_ptr<double> EshelbianPlasticity::OpBrokenAnalyticalTractionBc::piolaScalePtr

protected

Examples

EshelbianOperators.cpp.

Definition at line 573 of file EshelbianPlasticity.hpp.

scalingMethodsMap

std::map<std::string, boost::shared_ptr<ScalingMethod> > EshelbianPlasticity::OpBrokenAnalyticalTractionBc::scalingMethodsMap

protected

Definition at line 574 of file EshelbianPlasticity.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianPlasticity.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp