EshelbianPlasticity::dJ_dTractionImpl Struct Reference

Inheritance diagram for EshelbianPlasticity::dJ_dTractionImpl:

Collaboration diagram for EshelbianPlasticity::dJ_dTractionImpl:

Public Member Functions
MoFEMErrorCode	iNtegrate (EntData &data)
Public Member Functions inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
	OpAssembleBrokenFaceTopologicalDerivativeImplBase (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< TopologicalData > topo_ptr, boost::shared_ptr< double > J_ptr, SmartPetscObj< Vec > assemble_vec, Tag topo_tag, boost::shared_ptr< Range > ents_ptr=nullptr)
MoFEMErrorCode	aSsemble (EntData &data) override

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
using	OP = typename FormsIntegrators< FaceUserDataOperator >::Assembly< A >::OpBrokenBase
Protected Attributes inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
boost::shared_ptr< double >	JPtr
SmartPetscObj< Vec >	assembleVec
Tag	topoTag
boost::shared_ptr< TopologicalData >	topoData

Detailed Description

Definition at line 164 of file EshelbianTopologicalDerivativeOperators.cpp.

Member Function Documentation

iNtegrate()

MoFEMErrorCode EshelbianPlasticity::dJ_dTractionImpl::iNtegrate

(

EntData &

data

)

Definition at line 639 of file EshelbianTopologicalDerivativeOperators.cpp.

                                                        {
  MoFEMFunctionBegin;
 
#ifndef NDEBUG
  if (!topoData)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Topological data pointer is null");
#endif // NDEBUG
 
  const int nb_dofs = data.getIndices().size();
  if (!nb_dofs)
    MoFEMFunctionReturnHot(0);
 
  const int nb_integration_pts = OP::getGaussPts().size2();
 
#ifndef NDEBUG
  const int traction_rows = topoData->objDDisplacementAtPts.size1();
  const int traction_cols = topoData->objDDisplacementAtPts.size2();
  if (traction_cols != nb_integration_pts)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "objDDisplacementAtPts columns (%d) != nb integration points (%d)",
            traction_cols, nb_integration_pts);
 
  if (traction_rows != SPACE_DIM) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "objDDisplacementAtPts row size is %d, expected %d", traction_rows,
            SPACE_DIM);
  }
#endif // NDEBUG
 
  auto t_normal = OP::getFTensor1NormalsAtGaussPts();
  auto t_w = OP::getFTensor0IntegrationWeight();
  const int nb_base_functions = data.getN().size2() / SPACE_DIM;
  auto t_row_base_fun = data.getFTensor1N<SPACE_DIM>();
  auto t_obj_dtraction =
      getFTensor1FromMat<SPACE_DIM>(topoData->objDDisplacementAtPts);
 
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
 
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    auto t_nf = getFTensor1FromPtr<SPACE_DIM>(&*OP::locF.begin());
    int bb = 0;
    for (; bb != nb_dofs / SPACE_DIM; ++bb) {
      t_nf(i) +=
          t_w * (t_row_base_fun(j) * t_normal(j)) * t_obj_dtraction(i) * 0.5;
      ++t_nf;
      ++t_row_base_fun;
    }
    for (; bb != nb_base_functions; ++bb)
      ++t_row_base_fun;
 
    ++t_w;
    ++t_normal;
    ++t_obj_dtraction;
  }
 
  MoFEMFunctionReturn(0);
}

---

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

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