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

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

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

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

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

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

Detailed Description

template<int STRIDEMATD>
struct EshelbianPlasticity::HMHHencky::OpCalculateHelmholtzFreeEnergy< 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 185 of file HMHHencky.cpp.

Constructor & Destructor Documentation

OpCalculateHelmholtzFreeEnergy()

template<int STRIDEMATD>

EshelbianPlasticity::HMHHencky::OpCalculateHelmholtzFreeEnergy< STRIDEMATD >::OpCalculateHelmholtzFreeEnergy	(	boost::shared_ptr< DataAtIntegrationPts >	data_ptr,
		boost::shared_ptr< double >	total_helmholtz_free_energy_ptr
	)

Definition at line 1033 of file HMHHencky.cpp.

    : VolUserDataOperator(NOSPACE, OPSPACE), dataAtPts(data_ptr),
      totalHelmholtzFreeEnergyPtr(total_helmholtz_free_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::OpCalculateHelmholtzFreeEnergy< STRIDEMATD >::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 1048 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 (totalHelmholtzFreeEnergyPtr) {
    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;
    }
    *totalHelmholtzFreeEnergyPtr += getMeasure() * loc_energy;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

template<int STRIDEMATD>

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

private

Definition at line 193 of file HMHHencky.cpp.

totalHelmholtzFreeEnergyPtr

template<int STRIDEMATD>

boost::shared_ptr<double> EshelbianPlasticity::HMHHencky::OpCalculateHelmholtzFreeEnergy< STRIDEMATD >::totalHelmholtzFreeEnergyPtr

private

Definition at line 194 of file HMHHencky.cpp.

---

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

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