EshelbianPlasticity::HMHNeohookean::OpResidual Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHNeohookean::OpResidual:

Collaboration diagram for EshelbianPlasticity::HMHNeohookean::OpResidual:

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

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

Detailed Description

Definition at line 266 of file HMHNeohookean.cpp.

Constructor & Destructor Documentation

OpResidual()

EshelbianPlasticity::HMHNeohookean::OpResidual::OpResidual ( const std::string &  field, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, const double  alpha_u  )

inline

Definition at line 267 of file HMHNeohookean.cpp.

        : FormBase(field, field, FormBase::OPROW),
          dataAtPts(std::move(data_ptr)), alphaU(alpha_u) {
      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::OpResidual::iNtegrate ( EntData &  row_data )

inlineoverride

Definition at line 279 of file HMHNeohookean.cpp.

                                                         {
      MoFEMFunctionBegin;
      FTENSOR_INDEXES(3, i, j);
      FTENSOR_INDEX(size_symm, L);
      const auto t_packed_basis = FTensor::SymmLTensor<double, 3>();
      const auto material_ptr = getMaterial(dataAtPts);
      const auto parameters = material_ptr->getParameters(getFEEntityHandle());
      const double alpha_grad_u = material_ptr->alphaGradU;
      auto t_h = dataAtPts->getFTensorLogStretch(nbIntegrationPts);
      auto t_physical = dataAtPts->getFTensorAdjointPdU(nbIntegrationPts);
      // Inactive viscosity terms do not introduce rate-data dependencies.
      std::optional<decltype(dataAtPts->getFTensorLogStretchDot(
          nbIntegrationPts))>
          t_dot_h;
      std::optional<decltype(dataAtPts->getFTensorGradLogStretchDot(
          nbIntegrationPts))>
          t_grad_dot_h;
      if (alphaU != 0.)
        t_dot_h.emplace(dataAtPts->getFTensorLogStretchDot(nbIntegrationPts));
      if (alpha_grad_u != 0.)
        t_grad_dot_h.emplace(
            dataAtPts->getFTensorGradLogStretchDot(nbIntegrationPts));
      auto t_row = row_data.getFTensor0N();
      auto t_grad_row = row_data.getFTensor1DiffN<3>();
      auto t_weight = getFTensor0IntegrationWeight();
      for (int gg = 0; gg != nbIntegrationPts; ++gg) {
        Material::SymmetricTensor t_log_stretch;
        t_log_stretch(i, j) = t_h(i, j);
        PackedStress t_stress;
        double energy;
        CHKERR evaluateDirect(parameters, t_log_stretch, energy, &t_stress);
        const double alpha = getMeasure() * t_weight;
        PackedStress t_residual;
        t_residual(L) = t_stress(L) - t_physical(L);
        if (t_dot_h)
          t_residual(L) +=
              alphaU * (t_packed_basis(i, j, L) * (*t_dot_h)(i, j));
        FTensor::Tensor2<double, size_symm, 3> t_gradient;
        t_gradient(L, i) = 0.;
        if (t_grad_dot_h)
          t_gradient(L, i) = alpha_grad_u * (*t_grad_dot_h)(L, i);
        auto t_nf = getNf<size_symm>();
        int rr = 0;
        for (; rr != nbRows / size_symm; ++rr) {
          t_nf(L) += alpha *
                     (t_row * t_residual(L) + t_grad_row(i) * t_gradient(L, i));
          ++t_nf;
          ++t_row;
          ++t_grad_row;
        }
        for (; rr != nbRowBaseFunctions; ++rr) {
          ++t_row;
          ++t_grad_row;
        }
        ++t_weight;
        ++t_h;
        ++t_physical;
        if (t_dot_h)
          ++*t_dot_h;
        if (t_grad_dot_h)
          ++*t_grad_dot_h;
      }
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

alphaU

const double EshelbianPlasticity::HMHNeohookean::OpResidual::alphaU

private

Definition at line 346 of file HMHNeohookean.cpp.

dataAtPts

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

private

Definition at line 345 of file HMHNeohookean.cpp.

---

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

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