EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD > Struct Template Reference

Calculate energy density for Hencky material model. More...

Inheritance diagram for EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >:

Collaboration diagram for EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >:

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

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
boost::shared_ptr< double >	totalEnergyPtr

Detailed Description

template<int STRIDEMATD>
struct EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >

Calculate energy density for Hencky material model.

\[ \Psi(\log{\mathbf{U}}) = \frac{1}{2} U_{IJ} D_{IJKL} U_{KL} = \frac{1}{2} U_{IJ} T_{IJ} \]

where \(T_{IJ} = D_{IJKL} U_{KL}\) is a a Hencky stress.

Definition at line 218 of file HMHHencky.cpp.

Constructor & Destructor Documentation

OpCalculateEnergy()

template<int STRIDEMATD>

EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >::OpCalculateEnergy	(	boost::shared_ptr< DataAtIntegrationPts >	data_ptr,
		boost::shared_ptr< double >	total_energy_ptr
	)

Definition at line 1369 of file HMHHencky.cpp.

    : VolUserDataOperator(NOSPACE, OPSPACE), dataAtPts(data_ptr),
      totalEnergyPtr(total_energy_ptr) {
 
  if (!dataAtPts) {
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "dataAtPts is not allocated. Please set it before "
                      "using this operator.");
  }
}

Member Function Documentation

doWork()

template<int STRIDEMATD>

MoFEMErrorCode EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 1383 of file HMHHencky.cpp.

                                                                               {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, k);
  FTENSOR_INDEX(SPACE_DIM, l);
 
  int nb_integration_pts = getGaussPts().size2();
  auto t_log_u = dataAtPts->getFTensorLogStretchTotal(nb_integration_pts);
 
#ifndef NDEBUG
  auto &mat_d = dataAtPts->matD;
  if (mat_d.size2() != size_symm * size_symm) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "wrong matD size, number of columns should be %d but is %zu",
            size_symm * size_symm, mat_d.size2());
  }
  if constexpr (STRIDEMATD != 0) {
    if (mat_d.size1() != nb_integration_pts) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "wrong matD size, number of rows should be %d but is %zu",
              nb_integration_pts, mat_d.size1());
    }
  }
#endif
 
  auto t_D = getFTensor4DdgFromMat<3, 3, STRIDEMATD>(dataAtPts->matD);
 
  dataAtPts->energyAtPts.resize(nb_integration_pts, false);
  auto t_energy = getFTensor0FromVec(dataAtPts->energyAtPts);
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
 
    t_energy = 0.5 * (t_log_u(i, j) * (t_D(i, j, k, l) * t_log_u(k, l)));
 
    ++t_D;
    ++t_log_u;
    ++t_energy;
  }
 
  if (totalEnergyPtr) {
    auto t_w = getFTensor0IntegrationWeight();
    auto t_energy = getFTensor0FromVec(dataAtPts->energyAtPts);
    double loc_energy = 0;
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      loc_energy += t_energy * t_w;
      ++t_w;
      ++t_energy;
    }
    *totalEnergyPtr += getMeasure() * loc_energy;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

template<int STRIDEMATD>

boost::shared_ptr<DataAtIntegrationPts> EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >::dataAtPts

private

Definition at line 225 of file HMHHencky.cpp.

totalEnergyPtr

template<int STRIDEMATD>

boost::shared_ptr<double> EshelbianPlasticity::HMHHencky::OpCalculateEnergy< STRIDEMATD >::totalEnergyPtr

private

Definition at line 226 of file HMHHencky.cpp.

---

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

• users_modules/eshelbian_plasticity/src/impl/HMHHencky.cpp