OpFaceMaterialForce Struct Reference

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

Inheritance diagram for OpFaceMaterialForce:

Collaboration diagram for OpFaceMaterialForce:

Public Types
using	OP = FaceElementForcesAndSourcesCore::UserDataOperator

Public Member Functions
	OpFaceMaterialForce (boost::shared_ptr< DataAtIntegrationPts > data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts

Detailed Description

Definition at line 1002 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpFaceMaterialForce::OP = FaceElementForcesAndSourcesCore::UserDataOperator

Definition at line 1006 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpFaceMaterialForce()

OpFaceMaterialForce::OpFaceMaterialForce ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr )

inline

Definition at line 1008 of file EshelbianOperators.hpp.

      : OP(NOSPACE, OPSPACE), dataAtPts(data_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode OpFaceMaterialForce::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 4496 of file EshelbianOperators.cpp.

                                                          {
  MoFEMFunctionBegin;
 
#ifndef NDEBUG
  auto fe_mi_ptr = getFEMethod()->numeredEntFiniteElementPtr;
  auto pstatus = fe_mi_ptr->getPStatus();
  if (pstatus) {
    auto owner = fe_mi_ptr->getOwnerProc();
    MOFEM_LOG("SELF", Sev::noisy)
        << "OpFaceMaterialForce: owner proc is not 0, owner proc: " << owner
        << " " << getPtrFE()->mField.get_comm_rank();
  }
#endif // NDEBUG
 
  FTENSOR_INDEX(SPACE_DIM, I);
 
  FTensor::Tensor1<double, SPACE_DIM> t_face_T;
  t_face_T(I) = 0.;
  double face_pressure = 0.;
  auto t_T = dataAtPts->getFTensorFaceMaterialForce(
      getGaussPts().size2()); //< face material force
  auto t_p =
      getFTensor0FromVec(dataAtPts->normalPressureAtPts); //< normal pressure
  auto t_w = getFTensor0IntegrationWeight();
  for (auto gg = 0; gg != getGaussPts().size2(); ++gg) {
    t_face_T(I) += t_w * t_T(I);
    face_pressure += t_w * t_p;
    ++t_w;
    ++t_T;
    ++t_p;
  }
  t_face_T(I) *= getMeasure();
  face_pressure *= getMeasure();
 
  auto get_tag = [&](auto name, auto dim) {
    auto &moab = getPtrFE()->mField.get_moab();
    Tag tag;
    double def_val[] = {0., 0., 0.};
    CHK_MOAB_THROW(moab.tag_get_handle(name, dim, MB_TYPE_DOUBLE, tag,
                                       MB_TAG_CREAT | MB_TAG_SPARSE, def_val),
                   "create tag");
    return tag;
  };
 
  auto set_tag = [&](auto &&tag, auto ptr) {
    auto &moab = getPtrFE()->mField.get_moab();
    auto face = getPtrFE()->getFEEntityHandle();
    CHK_MOAB_THROW(moab.tag_set_data(tag, &face, 1, ptr), "set tag");
  };
 
  set_tag(get_tag("MaterialForce", 3), &t_face_T(0));
  set_tag(get_tag("FacePressure", 1), &face_pressure);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpFaceMaterialForce::dataAtPts

private

data at integration pts

Definition at line 1015 of file EshelbianOperators.hpp.

---

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

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