EshelbianPlasticity::HMHNeohookean::OpTangent Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHNeohookean::OpTangent:

Collaboration diagram for EshelbianPlasticity::HMHNeohookean::OpTangent:

Public Member Functions
	OpTangent (const std::string &row_field, const std::string &col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u)
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data) override

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
const double	alphaU

Detailed Description

Definition at line 356 of file HMHNeohookean.cpp.

Constructor & Destructor Documentation

OpTangent()

EshelbianPlasticity::HMHNeohookean::OpTangent::OpTangent ( const std::string &  row_field, const std::string &  col_field, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, const double  alpha_u  )

inline

Definition at line 357 of file HMHNeohookean.cpp.

        : FormBase(row_field, col_field, FormBase::OPROWCOL),
          dataAtPts(std::move(data_ptr)), alphaU(alpha_u) {
      sYmm = false;
      CHK_THROW_MESSAGE(checkDirectConfiguration(),
                        "Unsupported direct Neo-Hookean configuration");
      if (!std::isfinite(alphaU))
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                          "Neo-Hookean stretch viscosity must be finite");
    }

Member Function Documentation

iNtegrate()

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::OpTangent::iNtegrate ( EntData &  row_data, EntData &  col_data  )

inlineoverride

Definition at line 370 of file HMHNeohookean.cpp.

                                                                            {
      MoFEMFunctionBegin;
      FTENSOR_INDEXES(3, i, j, k, l);
      FTENSOR_INDEXES(size_symm, L, J);
      const auto material_ptr = getMaterial(dataAtPts);
      const auto parameters = material_ptr->getParameters(getFEEntityHandle());
      const auto t_packed_basis = FTensor::SymmLTensor<double, 3>();
      constexpr auto t_identity = FTensor::Kronecker_Delta<int>();
      PackedTangent t_viscous_metric;
      t_viscous_metric(L, J) =
          t_packed_basis(i, j, L) * t_packed_basis(i, j, J);
      const double ts_a = getTSa();
      const double alpha_grad_u = material_ptr->alphaGradU;
      auto t_h = dataAtPts->getFTensorLogStretch(nbIntegrationPts);
      auto t_pullback = dataAtPts->getFTensorAdjointPdstretch(nbIntegrationPts);
      auto t_eigenvalues = dataAtPts->getFTensorEigenVals(nbIntegrationPts);
      auto t_eigenvectors = dataAtPts->getFTensorEigenVecs(nbIntegrationPts);
      auto t_row = row_data.getFTensor0N();
      auto t_grad_row = row_data.getFTensor1DiffN<3>();
      auto t_weight = getFTensor0IntegrationWeight();
#ifndef NDEBUG
      if (dataAtPts->nbUniq.size() != nbIntegrationPts)
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Neo-Hookean tangent requires current geometric eigendata");
#endif
      for (int gg = 0; gg != nbIntegrationPts; ++gg) {
        Material::SymmetricTensor t_log_stretch;
        t_log_stretch(i, j) = t_h(i, j);
        PackedTangent t_tangent;
        double energy;
        CHKERR evaluateDirect(parameters, t_log_stretch, energy, nullptr,
                              &t_tangent);
        Material::SymmetricTensor t_symmetric_pullback;
        t_symmetric_pullback(i, j) =
            (t_pullback(i, j) || t_pullback(j, i)) / 2.;
        const auto t_work_curvature = EigenMatrix::getDiffDiffMat(
            t_eigenvalues, t_eigenvectors, EshelbianCore::f, EshelbianCore::d_f,
            EshelbianCore::dd_f, t_symmetric_pullback, dataAtPts->nbUniq[gg]);
        t_tangent(L, J) -=
            t_packed_basis(i, j, L) *
            (t_work_curvature(i, j, k, l) * t_packed_basis(k, l, J));
        t_tangent(L, J) += (alphaU * ts_a) * t_viscous_metric(L, J);
        const double alpha = getMeasure() * t_weight;
        int rr = 0;
        for (; rr != nbRows / size_symm; ++rr) {
          auto t_col = col_data.getFTensor0N(gg, 0);
          auto t_grad_col = col_data.getFTensor1DiffN<3>(gg, 0);
          auto t_matrix = getLocMat<size_symm>(size_symm * rr);
          for (int cc = 0; cc != nbCols / size_symm; ++cc) {
            t_matrix(L, J) += alpha * t_row * t_col * t_tangent(L, J);
            const double gradient_weight =
                alpha * alpha_grad_u * ts_a * (t_grad_row(i) * t_grad_col(i));
            t_matrix(L, J) += gradient_weight * t_identity(L, J);
            ++t_matrix;
            ++t_col;
            ++t_grad_col;
          }
          ++t_row;
          ++t_grad_row;
        }
        for (; rr != nbRowBaseFunctions; ++rr) {
          ++t_row;
          ++t_grad_row;
        }
        ++t_weight;
        ++t_h;
        ++t_pullback;
        ++t_eigenvalues;
        ++t_eigenvectors;
      }
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

alphaU

const double EshelbianPlasticity::HMHNeohookean::OpTangent::alphaU

private

Definition at line 445 of file HMHNeohookean.cpp.

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> EshelbianPlasticity::HMHNeohookean::OpTangent::dataAtPts

private

Definition at line 444 of file HMHNeohookean.cpp.

---

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

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