OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS > Struct Template Reference

#include <users_modules/multifield_plasticity/src/ElasticOperators.hpp>

Inheritance diagram for OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >:

Collaboration diagram for OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >:

Public Member Functions
	OpInternalForceRhs (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr)
	[Calculate stress] More...
MoFEMErrorCode	iNtegrate (EntData &data)
	[Internal force] More...
MoFEMErrorCode	aSsemble (EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

template<bool LOGSTRAIN, bool REACTIONS>
struct OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >

Definition at line 889 of file ElasticOperators.hpp.

Constructor & Destructor Documentation

OpInternalForceRhs()

template<bool LOGSTRAIN, bool REACTIONS>

OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >::OpInternalForceRhs	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

[Calculate stress]

Definition at line 392 of file ElasticOperators.cpp.

    : DomainEleOpAssembly(field_name, field_name, DomainEleOp::OPROW),
      commonDataPtr(common_data_ptr) {}

Member Function Documentation

aSsemble()

template<bool LOGSTRAIN, bool REACTIONS>

MoFEMErrorCode OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >::aSsemble

(

EntData &

data

)

Definition at line 449 of file ElasticOperators.cpp.

                                                                {
  MoFEMFunctionBeginHot;
  if (REACTIONS) {
    auto react_ents = commonDataPtr->reactionEnts;
    EntityHandle ent = getFEEntityHandle();
 
    auto &vec = commonDataPtr->reactionsVec;
    double *react_ptr;
    CHKERR VecGetArray(vec, &react_ptr);
 
    int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
    int ii = 0;
    for (auto &dof : data.getFieldDofs()) {
      auto ent = dof->getEnt();
      double tag_val[rank];
      auto idx = dof->getDofCoeffIdx();
      auto &m_field = getPtrFE()->mField;
      CHKERR m_field.get_moab().tag_get_data(commonDataPtr->reactionTag, &ent,
                                             1, tag_val);
      if (dof->getEntType() == MBVERTEX) {
        const double react = this->locF(ii) / (*cache).young_modulus_inv;
        tag_val[idx] += react;
        if (react_ents.find(ent) != react_ents.end())
          react_ptr[idx] += react;
      }
      ii++;
      CHKERR m_field.get_moab().tag_set_data(commonDataPtr->reactionTag, &ent,
                                             1, &tag_val);
    }
    CHKERR VecRestoreArray(vec, &react_ptr);
  } else
    CHKERR VecSetValues<EssentialBcStorage>(
        this->getKSPf(), data, &*this->locF.data().begin(), ADD_VALUES);
 
  MoFEMFunctionReturnHot(0);
}

iNtegrate()

template<bool LOGSTRAIN, bool REACTIONS>

MoFEMErrorCode OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >::iNtegrate

(

EntData &

data

)

[Internal force]

Definition at line 400 of file ElasticOperators.cpp.

                                                                 {
  MoFEMFunctionBegin;
 
  const size_t nb_dofs = data.getIndices().size();
  if (nb_dofs) {
 
    const size_t nb_base_functions = data.getN().size2();
    if (3 * nb_base_functions < nb_dofs)
      SETERRQ(
          PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
          "Number of DOFs is larger than number of base functions on entity");
 
    const size_t nb_gauss_pts = data.getN().size1();
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_stress =
        getFTensor2SymmetricFromMat<3>(*(commonDataPtr->mStressPtr));
    auto t_piola = getFTensor2FromMat<3, 3>(*(commonDataPtr->mPiolaStressPtr));
    auto t_diff_base = data.getFTensor1DiffN<3>();
 
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
      double alpha = getMeasure() * t_w;
      auto t_nf = getFTensor1FromArray<3, 3>(locF);
      auto out_stress = to_non_symm(t_stress);
 
      if (LOGSTRAIN) {
        out_stress(i, j) = t_piola(i, j);
        ++t_piola;
      }
 
      size_t bb = 0;
      for (; bb != nb_dofs / 3; ++bb) {
        t_nf(i) += alpha * t_diff_base(j) * out_stress(i, j);
        ++t_diff_base;
        ++t_nf;
      }
      for (; bb < nb_base_functions; ++bb)
        ++t_diff_base;
 
      ++t_stress;
      ++t_w;
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

template<bool LOGSTRAIN, bool REACTIONS>

boost::shared_ptr<CommonData> OpElasticTools::OpInternalForceRhs< LOGSTRAIN, REACTIONS >::commonDataPtr

private

Definition at line 896 of file ElasticOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/ElasticOperators.hpp
• users_modules/multifield_plasticity/src/impl/ElasticOperators.cpp