EshelbianPlasticity::OpAnalyticalDispBc_dX Struct Reference

Inheritance diagram for EshelbianPlasticity::OpAnalyticalDispBc_dX:

Collaboration diagram for EshelbianPlasticity::OpAnalyticalDispBc_dX:

Public Member Functions
	OpAnalyticalDispBc_dX (const std::string &field_name, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< AnalyticalDisplacementBcVec > bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, boost::shared_ptr< TopologicalData > topo_ptr, SmartPetscObj< Vec > vec, boost::shared_ptr< double > J_ptr=nullptr, Tag tag=Tag())
MoFEMErrorCode	integrate (int side, EntityType type, EntData &data) override
Public Member Functions inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
	OpAssembleTopologicalObjectiveDerivativeImplBase (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< TopologicalData > topo_ptr, boost::shared_ptr< double > J_ptr, SmartPetscObj< Vec > assemble_vec, Tag topo_tag)
MoFEMErrorCode	assemble (int side, EntityType type, EntData &data) override
Public Member Functions inherited from OpAssembleBasic< FaceUserDataOperator >
	OpAssembleBasic (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)
	OpAssembleBasic (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry, ScaleOff scale_off=[]() { return 1;})
	OpAssembleBasic (const FieldSpace space)
virtual MoFEMErrorCode	integrate (EntData &data)
virtual MoFEMErrorCode	integrate (int row_side, EntityType row_type, EntData &data)
virtual MoFEMErrorCode	integrate (EntData &row_data, EntData &col_data)
virtual MoFEMErrorCode	assemble (EntData &data)
virtual MoFEMErrorCode	assemble (int row_side, EntityType row_type, EntData &data)
virtual MoFEMErrorCode	assemble (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

Private Attributes
boost::shared_ptr< std::vector< BrokenBaseSideData > >	brokenSideDataPtr
boost::shared_ptr< AnalyticalDisplacementBcVec >	bcDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
using	OP = OpAssembleFace
Public Types inherited from OpAssembleBasic< FaceUserDataOperator >
using	ScaleOff = boost::function< double()>
Public Attributes inherited from OpAssembleBasic< FaceUserDataOperator >
const bool	assembleSymmetry
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts
VectorDouble	nF
	local right hand side vector
MatrixDouble	K
	local tangent matrix
MatrixDouble	transposeK
ScaleOff	scaleOff
Protected Attributes inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
double	locJ
boost::shared_ptr< double >	JPtr
SmartPetscObj< Vec >	assembleVec
Tag	topoTag
boost::shared_ptr< TopologicalData >	topoData

Detailed Description

Definition at line 256 of file EshelbianTopologicalDerivativeOperators.cpp.

Constructor & Destructor Documentation

OpAnalyticalDispBc_dX()

EshelbianPlasticity::OpAnalyticalDispBc_dX::OpAnalyticalDispBc_dX ( const std::string &  field_name, boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_side_data_ptr, boost::shared_ptr< AnalyticalDisplacementBcVec >  bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > >  smv, boost::shared_ptr< TopologicalData >  topo_ptr, SmartPetscObj< Vec >  vec, boost::shared_ptr< double >  J_ptr = nullptr, Tag  tag = Tag()  )

inline

Definition at line 257 of file EshelbianTopologicalDerivativeOperators.cpp.

      : OpAssembleFaceTopologicalDerivativeImpl(field_name, nullptr, topo_ptr,
                                                J_ptr, vec, tag),
        brokenSideDataPtr(broken_side_data_ptr), bcDispPtr(bc_disp_ptr),
        scalingMethodsMap(smv) {}

Member Function Documentation

integrate()

MoFEMErrorCode EshelbianPlasticity::OpAnalyticalDispBc_dX::integrate ( int  side, EntityType  type, EntData &  data  )

override

Definition at line 1285 of file EshelbianTopologicalDerivativeOperators.cpp.

                                                               {
  MoFEMFunctionBegin;
  locJ = 0;
 
#ifndef NDEBUG
  if (!brokenSideDataPtr)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Broken side data pointer is null");
#endif // NDEBUG
 
  const int nb_dofs = data.getIndices().size();
  if (!nb_dofs)
    MoFEMFunctionReturnHot(0);
 
  const int nb_integration_pts = getGaussPts().size2();
  const int nb_base_functions = data.getN().size2();
 
  if (data.getDiffN().size1() != nb_integration_pts)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Differential of base functions should have the same number of "
            "integration points as the data");
  if (data.getDiffN().size2() != nb_base_functions * 2)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Differential of base functions should have the same number of "
            "base functions as the data");
 
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, k);
 
  const EntityHandle fe_ent = getFEEntityHandle();
  for (auto &bc : *bcDispPtr) {
    if (bc.faces.find(fe_ent) == bc.faces.end())
      continue;
 
    auto v_analytical_expr =
        getTopologicalAnalyticalExpr(this, bc.blockName);
 
    for (auto &bd : *brokenSideDataPtr) {
      auto t_w = getFTensor0IntegrationWeight();
      auto t_tangent1 = getFTensor1Tangent1AtGaussPts();
      auto t_tangent2 = getFTensor1Tangent2AtGaussPts();
      auto t_var_flux =
          getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(bd.getVarFlux());
      auto t_diff_base = data.getFTensor1DiffN<2>();
      auto t_bc_disp =
          getFTensor1FromMat<SPACE_DIM, -1, DL>(v_analytical_expr);
 
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        const double a = 0.5 * bd.getSense() * t_w;
 
        FTensor::Tensor1<double, SPACE_DIM> t_normal;
        t_normal(j) =
            (FTensor::levi_civita(i, j, k) * t_tangent1(k)) * t_tangent2(i);
 
        locJ += a * t_bc_disp(i) * (t_var_flux(i, j) * t_normal(j));
 
        auto t_nf = getFTensor1FromArray<SPACE_DIM, SPACE_DIM>(nF);
        int bb = 0;
        for (; bb != nb_dofs / SPACE_DIM; ++bb) {
          // The H(div) contravariant Piola trace P . n dA is invariant with
          // respect to material-position perturbations. The normal variation is
          // cancelled by the Piola variation of the flux, so this BC contributes
          // to locJ but not to the exact dX vector.
          ++t_nf;
          ++t_diff_base;
        }
        for (; bb != nb_base_functions; ++bb)
          ++t_diff_base;
 
        ++t_w;
        ++t_tangent1;
        ++t_tangent2;
        ++t_var_flux;
        ++t_bc_disp;
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcDispPtr

boost::shared_ptr<AnalyticalDisplacementBcVec> EshelbianPlasticity::OpAnalyticalDispBc_dX::bcDispPtr

private

Definition at line 273 of file EshelbianTopologicalDerivativeOperators.cpp.

brokenSideDataPtr

boost::shared_ptr<std::vector<BrokenBaseSideData> > EshelbianPlasticity::OpAnalyticalDispBc_dX::brokenSideDataPtr

private

Definition at line 272 of file EshelbianTopologicalDerivativeOperators.cpp.

scalingMethodsMap

std::map<std::string, boost::shared_ptr<ScalingMethod> > EshelbianPlasticity::OpAnalyticalDispBc_dX::scalingMethodsMap

private

Definition at line 274 of file EshelbianTopologicalDerivativeOperators.cpp.

---

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

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