OpReleaseEnergy Struct Reference

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

Inheritance diagram for OpReleaseEnergy:

Collaboration diagram for OpReleaseEnergy:

Public Types
using	OP = VolumeElementForcesAndSourcesCore::UserDataOperator

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

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts More...
boost::shared_ptr< double >	releaseEnergyPtr

Detailed Description

Examples

EshelbianPlasticity.cpp.

Definition at line 577 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpReleaseEnergy::OP = VolumeElementForcesAndSourcesCore::UserDataOperator

Definition at line 580 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpReleaseEnergy()

OpReleaseEnergy::OpReleaseEnergy ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr, boost::shared_ptr< double >  release_energy_ptr  )

inline

Definition at line 582 of file EshelbianOperators.hpp.

      : OP(NOSPACE, OPSPACE), dataAtPts(data_ptr),
        releaseEnergyPtr(release_energy_ptr) {
    if (!releaseEnergyPtr)
      CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                        "releaseEnergyPtr is nullptr");
  }

Member Function Documentation

doWork()

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

Examples

EshelbianOperators.cpp.

Definition at line 2436 of file EshelbianOperators.cpp.

                                                      {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, I);
  FTENSOR_INDEX(SPACE_DIM, J);
  FTENSOR_INDEX(SPACE_DIM, K);
 
  // sym size indices
  FTENSOR_INDEX(size_symm, L);
  auto t_L = symm_L_tensor();
 
  const auto nb_integration_pts = getGaussPts().size2();
 
  auto t_w = getFTensor0IntegrationWeight();
  auto t_h = getFTensor2FromMat<3, 3>(dataAtPts->hAtPts);
  auto t_P = getFTensor2FromMat<3, 3>(dataAtPts->approxPAtPts);
 
  auto t_var_omega = getFTensor1FromMat<3>(dataAtPts->varRotAxis);
  auto t_var_log_u = getFTensor2SymmetricFromMat<3>(dataAtPts->varLogStreach);
  auto t_var_P = getFTensor2FromMat<3, 3>(dataAtPts->varPiola);
 
  auto t_h_domega = getFTensor3FromMat<3, 3, 3>(dataAtPts->hdOmegaAtPts);
  auto t_h_dlog_u =
      getFTensor3FromMat<3, 3, size_symm>(dataAtPts->hdLogStretchAtPts);
 
  auto t_approx_P_adjoint_log_du =
      getFTensor1FromMat<size_symm>(dataAtPts->adjointPdUAtPts);
 
  double var_element_energy = 0.;
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
 
    FTensor::Tensor1<double, size_symm> t_var_L_u;
    t_var_L_u(L) = t_L(J, K, L) * t_var_log_u(J, K);
    auto var_energy = t_P(i, J) * (t_h_domega(i, J, j) * t_var_omega(j) +
                                   t_h_dlog_u(i, J, L) * t_var_L_u(L));
    var_element_energy += t_w * var_energy;
    auto var_complementary = t_var_P(i, J) * t_h(i, J);
    var_element_energy += t_w * var_complementary;
 
    ++t_w;
    ++t_h;
    ++t_P;
 
    ++t_var_omega;
    ++t_var_log_u;
    ++t_var_P;
 
    ++t_h_domega;
    ++t_h_dlog_u;
 
    ++t_approx_P_adjoint_log_du;
  }
 
  var_element_energy *= getMeasure();
 
  auto get_tag = [&]() {
    auto &mob = getPtrFE()->mField.get_moab();
    Tag tag;
    double def_val[] = {0.};
    CHK_MOAB_THROW(mob.tag_get_handle("ReleaseEnergy", 1, MB_TYPE_DOUBLE, tag,
                                      MB_TAG_CREAT | MB_TAG_SPARSE, def_val),
                   "create tag");
    return tag;
  };
 
  auto set_tag = [&](auto &&tag, auto &energy) {
    auto &mob = getPtrFE()->mField.get_moab();
    auto tet = getPtrFE()->getFEEntityHandle();
    CHK_MOAB_THROW(mob.tag_set_data(tag, &tet, 1, &energy), "set tag");
  };
 
  set_tag(get_tag(), var_element_energy);
 
  *releaseEnergyPtr += var_element_energy;
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpReleaseEnergy::dataAtPts

private

data at integration pts

Definition at line 595 of file EshelbianOperators.hpp.

releaseEnergyPtr

boost::shared_ptr<double> OpReleaseEnergy::releaseEnergyPtr

private

Definition at line 596 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