EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress:

Collaboration diagram for EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress:

Public Member Functions
	OpCalculateStretchFromStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< HMHHencky > hencky_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts
boost::shared_ptr< MatrixDouble >	strainPtr
boost::shared_ptr< MatrixDouble >	stressPtr
boost::shared_ptr< HMHHencky >	henckyPtr

Detailed Description

Definition at line 177 of file HMHHencky.cpp.

Constructor & Destructor Documentation

OpCalculateStretchFromStress()

EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress::OpCalculateStretchFromStress	(	boost::shared_ptr< DataAtIntegrationPts >	data_ptr,
		boost::shared_ptr< MatrixDouble >	strain_ptr,
		boost::shared_ptr< MatrixDouble >	stress_ptr,
		boost::shared_ptr< HMHHencky >	hencky_ptr )

Definition at line 995 of file HMHHencky.cpp.

    : VolUserDataOperator(NOSPACE, OPSPACE), dataAtPts(data_ptr),
      strainPtr(strain_ptr), stressPtr(stress_ptr), henckyPtr(hencky_ptr) {}

Member Function Documentation

doWork()

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

Definition at line 1003 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);
  FTENSOR_INDEX(SPACE_DIM, m);
  FTENSOR_INDEX(SPACE_DIM, n);
 
  auto nb_integration_pts = stressPtr->size2();
#ifndef NDEBUG
  if (nb_integration_pts != getGaussPts().size2()) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "inconsistent number of integration points");
  }
#endif // NDEBUG
 
  auto get_D = [&]() {
    MoFEMFunctionBegin;
    for (auto &b : henckyPtr->blockData) {
 
      if (b.blockEnts.find(getFEEntityHandle()) != b.blockEnts.end()) {
        dataAtPts->mu = b.shearModulusG;
        dataAtPts->lambda = b.bulkModulusK - 2 * b.shearModulusG / 3;
        CHKERR henckyPtr->getMatDPtr(
            dataAtPts->getMatDPtr(), dataAtPts->getMatAxiatorDPtr(),
            dataAtPts->getMatDeviatorDPtr(), b.bulkModulusK, b.shearModulusG);
        CHKERR henckyPtr->getInvMatDPtr(dataAtPts->getMatInvDPtr(),
                                        b.bulkModulusK, b.shearModulusG);
        MoFEMFunctionReturnHot(0);
      }
    }
 
    const auto E = henckyPtr->E;
    const auto nu = henckyPtr->nu;
 
    double bulk_modulus_K = E / (3 * (1 - 2 * nu));
    double shear_modulus_G = E / (2 * (1 + nu));
 
    dataAtPts->mu = shear_modulus_G;
    dataAtPts->lambda = bulk_modulus_K - 2 * shear_modulus_G / 3;
 
    CHKERR henckyPtr->getMatDPtr(
        dataAtPts->getMatDPtr(), dataAtPts->getMatAxiatorDPtr(),
        dataAtPts->getMatDeviatorDPtr(), bulk_modulus_K, shear_modulus_G);
    CHKERR henckyPtr->getInvMatDPtr(dataAtPts->getMatInvDPtr(), bulk_modulus_K,
                                    shear_modulus_G);
    MoFEMFunctionReturn(0);
  };
 
  CHKERR get_D();
 
  strainPtr->resize(size_symm, nb_integration_pts, false);
  auto t_strain = getFTensor2SymmetricFromMat<3>(*strainPtr);
  auto t_stress = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*stressPtr);
  auto t_inv_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
      &*dataAtPts->matInvD.data().begin());
#ifndef NDEBUG
  auto t_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
      &*dataAtPts->matD.data().begin());
#endif
 
  const double lambda = dataAtPts->lambda;
  const double mu = dataAtPts->mu;
  constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
 
  // note: add rotation, so we can extract rigid body motion, work then with
  // symmetric part.
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
    t_strain(i, j) = t_inv_D(i, j, k, l) * t_stress(k, l);
 
#ifndef NDEBUG
    FTensor::Tensor2_symmetric<double, 3> t_stress_symm_debug;
    t_stress_symm_debug(i, j) = (t_stress(i, j) || t_stress(j, i)) / 2;
    FTensor::Tensor2_symmetric<double, 3> t_stress_symm_debug_diff;
    t_stress_symm_debug_diff(i, j) =
        t_D(i, j, k, l) * t_strain(k, l) - t_stress_symm_debug(i, j);
    double nrm =
        t_stress_symm_debug_diff(i, j) * t_stress_symm_debug_diff(i, j);
    double nrm0 = t_stress_symm_debug(i, j) * t_stress_symm_debug(i, j) +
                  std::numeric_limits<double>::epsilon();
    constexpr double eps = 1e-10;
    if (std::fabs(std::sqrt(nrm / nrm0)) > eps) {
      MOFEM_LOG("SELF", Sev::error)
          << "Stress symmetry check failed: " << std::endl
          << t_stress_symm_debug_diff << std::endl
          << t_stress;
      CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                        "Norm is too big: " + std::to_string(nrm / nrm0));
    }
    ++t_D;
#endif
 
    ++t_strain;
    ++t_stress;
    ++t_inv_D;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

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

private

data at integration pts

Definition at line 187 of file HMHHencky.cpp.

henckyPtr

boost::shared_ptr<HMHHencky> EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress::henckyPtr

private

Definition at line 190 of file HMHHencky.cpp.

strainPtr

boost::shared_ptr<MatrixDouble> EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress::strainPtr

private

Definition at line 188 of file HMHHencky.cpp.

stressPtr

boost::shared_ptr<MatrixDouble> EshelbianPlasticity::HMHHencky::OpCalculateStretchFromStress::stressPtr

private

Definition at line 189 of file HMHHencky.cpp.

---

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

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