EshelbianPlasticity::HMHHencky::OpCalculateEnergy Struct Reference

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

Inheritance diagram for EshelbianPlasticity::HMHHencky::OpCalculateEnergy:

Collaboration diagram for EshelbianPlasticity::HMHHencky::OpCalculateEnergy:

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

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 160 of file HMHHencky.cpp.

Constructor & Destructor Documentation

OpCalculateEnergy()

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

Definition at line 932 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()

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

Definition at line 946 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 =
      getFTensor2SymmetricFromMat<3>(dataAtPts->logStretchTotalTensorAtPts);
 
 
  auto &mat_d = dataAtPts->matD;
  if (mat_d.size1() != size_symm || mat_d.size2() != size_symm) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "wrong matD size, should be 6 by 6 but is %d by %d", size_symm,
             size_symm);
  }
 
  auto t_D = getFTensor4DdgFromPtr<3, 3, 0>(mat_d.data().data());
 
  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_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

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

private

Definition at line 167 of file HMHHencky.cpp.

totalEnergyPtr

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

private

Definition at line 168 of file HMHHencky.cpp.

---

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

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