OpGetInternalStress< INTERP_ORDER > Struct Template Reference

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

Inheritance diagram for OpGetInternalStress< INTERP_ORDER >:

Collaboration diagram for OpGetInternalStress< INTERP_ORDER >:

Public Member Functions
	OpGetInternalStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, std::string tag_name="")
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
std::string	tagName

Detailed Description

template<int INTERP_ORDER>
struct OpGetInternalStress< INTERP_ORDER >

Examples

EshelbianPlasticity.cpp.

Definition at line 392 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpGetInternalStress()

template<int INTERP_ORDER>

OpGetInternalStress< INTERP_ORDER >::OpGetInternalStress ( boost::shared_ptr< DataAtIntegrationPts > data_ptr, std::string tag_name = "" )

inline

Definition at line 393 of file EshelbianOperators.hpp.

      : UserDataOperator(H1, OPLAST), dataAtPts(data_ptr), tagName(tag_name) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBVERTEX] = true;
  }

Member Function Documentation

doWork() [1/3]

template<int INTERP_ORDER>

MoFEMErrorCode OpGetInternalStress< INTERP_ORDER >::doWork	(	int	side,
		EntityType	type,
		EntData &	data )

doWork() [2/3]

MoFEMErrorCode OpGetInternalStress< 0 >::doWork	(	int	side,
		EntityType	type,
		EntData &	data )

Definition at line 1241 of file EshelbianOperators.cpp.

                                                             {
  MoFEMFunctionBegin;
 
  int nb_integration_pts = getGaussPts().size2();
 
  Tag tag;
  CHKERR getPtrFE() -> mField.get_moab().tag_get_handle(tagName.c_str(), tag);
  int tag_length;
  CHKERR getPtrFE() -> mField.get_moab().tag_get_length(tag, tag_length);
  if (tag_length != 9) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Number of internal stress components should be 9 but is %d",
             tag_length);
  }
 
  VectorDouble const_stress_vec(9);
  auto fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
  CHKERR getPtrFE() -> mField.get_moab().tag_get_data(
      tag, &fe_ent, 1, &*const_stress_vec.data().begin());
  auto t_const_stress = getFTensor1FromArray<9, 9>(const_stress_vec);
 
  dataAtPts->internalStressAtPts.resize(tag_length, nb_integration_pts, false);
  dataAtPts->internalStressAtPts.clear();
  auto t_internal_stress =
      getFTensor1FromMat<9>(dataAtPts->internalStressAtPts);
 
  FTensor::Index<'L', 9> L;
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    t_internal_stress(L) = t_const_stress(L);
    ++t_internal_stress;
  }
 
  MoFEMFunctionReturn(0);
}

doWork() [3/3]

MoFEMErrorCode OpGetInternalStress< 1 >::doWork	(	int	side,
		EntityType	type,
		EntData &	data )

Definition at line 1278 of file EshelbianOperators.cpp.

                                                                      {
  MoFEMFunctionBegin;
 
  int nb_integration_pts = getGaussPts().size2();
 
  Tag tag;
  CHKERR getPtrFE() -> mField.get_moab().tag_get_handle(tagName.c_str(), tag);
  int tag_length;
  CHKERR getPtrFE() -> mField.get_moab().tag_get_length(tag, tag_length);
  if (tag_length != 9) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Number of internal stress components should be 9 but is %d",
             tag_length);
  }
 
  auto fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
  const EntityHandle *vert_conn;
  int vert_num;
  CHKERR getPtrFE() -> mField.get_moab().get_connectivity(fe_ent, vert_conn,
                                                          vert_num, true);
  VectorDouble vert_data(vert_num * tag_length);
  CHKERR getPtrFE() -> mField.get_moab().tag_get_data(tag, vert_conn, vert_num,
                                                      &vert_data[0]);
 
  dataAtPts->internalStressAtPts.resize(tag_length, nb_integration_pts, false);
  dataAtPts->internalStressAtPts.clear();
  auto t_internal_stress =
      getFTensor1FromMat<9>(dataAtPts->internalStressAtPts);
 
  auto t_shape_n = data.getFTensor0N();
  int nb_shape_fn = data.getN(NOBASE).size2();
  FTensor::Index<'L', 9> L;
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    auto t_vert_data = getFTensor1FromArray<9, 9>(vert_data);
    for (int bb = 0; bb != nb_shape_fn; ++bb) {
      t_internal_stress(L) += t_vert_data(L) * t_shape_n;
      ++t_vert_data;
      ++t_shape_n;
    }
    ++t_internal_stress;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

template<int INTERP_ORDER>

boost::shared_ptr<DataAtIntegrationPts> OpGetInternalStress< INTERP_ORDER >::dataAtPts

private

Definition at line 402 of file EshelbianOperators.hpp.

tagName

template<int INTERP_ORDER>

std::string OpGetInternalStress< INTERP_ORDER >::tagName

private

Definition at line 403 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp