Inheritance diagram for EshelbianPlasticity::OpDispBc_dX:

Collaboration diagram for EshelbianPlasticity::OpDispBc_dX:

Public Member Functions
	OpDispBc_dX (const std::string &field_name, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< BcDispVec > 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< BcDispVec >	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 205 of file EshelbianTopologicalDerivativeOperators.cpp.

Constructor & Destructor Documentation

◆ OpDispBc_dX()

EshelbianPlasticity::OpDispBc_dX::OpDispBc_dX	(	const std::string &	field_name,
		boost::shared_ptr< std::vector< BrokenBaseSideData > >	broken_side_data_ptr,
		boost::shared_ptr< BcDispVec >	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 206 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::OpDispBc_dX::integrate	(	int	side,
		EntityType	type,
		EntData &	data
	)

override

Definition at line 1189 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");
 
  double time = getFEMethod()->ts_t;
  if (EshelbianCore::physicalTimeFlg)
    time = EshelbianCore::currentPhysicalTime;
 
  FTENSOR_INDEXES(SPACE_DIM, i, j, k, I);
 
  FTensor::Number<0> N0;
  FTensor::Number<1> N1;
 
  const EntityHandle fe_ent = getFEEntityHandle();
  for (auto &bc : *bcDispPtr) {
    if (bc.faces.find(fe_ent) == bc.faces.end())
      continue;
 
    double scale = 1;
    if (scalingMethodsMap.find(bc.blockName) != scalingMethodsMap.end()) {
      scale *= scalingMethodsMap.at(bc.blockName)->getScale(time);
    } else {
      MOFEM_LOG("SELF", Sev::warning)
          << "No scaling method found for " << bc.blockName;
    }
 
    FTensor::Tensor1<double, SPACE_DIM> t_bc_disp(bc.vals[0], bc.vals[1],
                                                  bc.vals[2]);
    t_bc_disp(i) *= scale;
 
    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>();
 
      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;
      }
    }
 
  }
 
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

◆ bcDispPtr

boost::shared_ptr<BcDispVec> EshelbianPlasticity::OpDispBc_dX::bcDispPtr

private

Definition at line 222 of file EshelbianTopologicalDerivativeOperators.cpp.

◆ brokenSideDataPtr

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

private

Definition at line 221 of file EshelbianTopologicalDerivativeOperators.cpp.

◆ scalingMethodsMap

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

private

Definition at line 223 of file EshelbianTopologicalDerivativeOperators.cpp.

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

users_modules/eshelbian_plasticity/src/impl/EshelbianTopologicalDerivativeOperators.cpp

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ OpDispBc_dX()

Member Function Documentation

◆ integrate()

Member Data Documentation

◆ bcDispPtr

◆ brokenSideDataPtr

◆ scalingMethodsMap