OpTauStabilizationDispRhsBc Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpTauStabilizationDispRhsBc:

Collaboration diagram for OpTauStabilizationDispRhsBc:

Public Types
using	OP = OpBrokenBaseTimesBrokenDisp
Public Types inherited from OpBrokenBaseTimesBrokenDisp
using	OP = OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
Public Types inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
using	BASE = OpBrokenBaseImpl< OpBaseTimesVectorFace >

Public Member Functions
	OpTauStabilizationDispRhsBc (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< BcDispVec > &bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, ScalarFun tau_coeff, boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from OpBrokenBaseTimesBrokenDisp
	OpBrokenBaseTimesBrokenDisp (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, ScalarFun beta_coeff=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &data)
Protected Member Functions inherited from OpBrokenBaseTimesBrokenDisp
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Protected Attributes
boost::shared_ptr< BcDispVec >	bcDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap
Protected Attributes inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
boost::weak_ptr< MatrixDouble >	fluxMatPtr

Detailed Description

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 1246 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpTauStabilizationDispRhsBc::OP = OpBrokenBaseTimesBrokenDisp

Definition at line 1248 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpTauStabilizationDispRhsBc()

OpTauStabilizationDispRhsBc::OpTauStabilizationDispRhsBc ( boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data, boost::shared_ptr< BcDispVec > &  bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > >  smv, ScalarFun  tau_coeff, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 1250 of file EshelbianOperators.hpp.

      : OpBrokenBaseTimesBrokenDisp(broken_base_side_data, tau_coeff, ents_ptr),
        bcDispPtr(bc_disp_ptr), scalingMethodsMap(smv) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpTauStabilizationDispRhsBc::iNtegrate ( EntData &  data )

protected

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1613 of file EshelbianOperators.cpp.

                                                                   {
  MoFEMFunctionBegin;
  // get entity of face
  EntityHandle fe_ent = OP::getFEEntityHandle();
  // iterate over all boundary data
  for (auto &bc : (*bcDispPtr)) {
    // check if finite element entity is part of boundary condition
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
      int nb_dofs = data.getIndices().size();
 
      int nb_integration_pts = OP::getGaussPts().size2();
      auto t_w = OP::getFTensor0IntegrationWeight();
      int nb_base_functions = data.getN().size2();
      auto t_row_base_fun = data.getFTensor0N();
#ifndef NDEBUG
      if (!this->sourceVec) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Source vector for OpTauStabilizationDispRhsBc is not set");
      }
      if (data.getN().size1() != nb_integration_pts) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Number of integration points in data should be %d but is %d",
                nb_integration_pts, (int)data.getN().size1());
      }
      if (nb_base_functions < nb_dofs / SPACE_DIM) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Number of base functions in data should be %d but is %d",
                nb_base_functions, (int)data.getN().size2() / SPACE_DIM);
      }
 
#endif
 
      auto t_disp_val =
          MatrixSizeHelper<GetFTensor1FromMatType<3, -1, DL>, DL>::get(
              *this->sourceVec, nb_integration_pts)();
 
      FTENSOR_INDEX(SPACE_DIM, i);
      FTENSOR_INDEX(SPACE_DIM, j);
 
      double scale = 1;
      if (scalingMethodsMap.find(bc.blockName) != scalingMethodsMap.end()) {
        if (EshelbianCore::physicalTimeFlg) {
          scale *= scalingMethodsMap.at(bc.blockName)
                       ->getScale(EshelbianCore::currentPhysicalTime);
        } else {
          scale *= scalingMethodsMap.at(bc.blockName)
                       ->getScale(OP::getFEMethod()->ts_t);
        }
      } else {
        MOFEM_LOG("SELF", Sev::warning)
            << "No scaling method found for " << bc.blockName;
      }
 
      // get bc data
      FTensor::Tensor1<double, 3> t_bc_disp(bc.vals[0], bc.vals[1], bc.vals[2]);
      t_bc_disp(i) *= scale;
 
      auto area = getMeasure();
      auto t_coords = OP::getFTensor1CoordsAtGaussPts();
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        auto tau_scale =
            area * t_w * OP::betaCoeff(t_coords(0), t_coords(1), t_coords(2));
        FTensor::Tensor1<double, SPACE_DIM> t_bc_residual;
        for (int ii = 0; ii != SPACE_DIM; ++ii)
          t_bc_residual(ii) =
              bc.flags[ii] ? t_disp_val(ii) - t_bc_disp(ii) : 0.;
        auto t_nf = getFTensor1FromPtr<3, 3>(OP::locF.data().data());
        int bb = 0;
        for (; bb != nb_dofs / SPACE_DIM; ++bb) {
          t_nf(i) +=
              (tau_scale * t_row_base_fun) * t_bc_residual(i);
          ++t_nf;
          ++t_row_base_fun;
        }
        for (; bb != nb_base_functions; ++bb)
          ++t_row_base_fun;
 
        ++t_w;
        ++t_coords;
        ++t_disp_val;
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcDispPtr

boost::shared_ptr<BcDispVec> OpTauStabilizationDispRhsBc::bcDispPtr

protected

Definition at line 1260 of file EshelbianOperators.hpp.

scalingMethodsMap

std::map<std::string, boost::shared_ptr<ScalingMethod> > OpTauStabilizationDispRhsBc::scalingMethodsMap

protected

Definition at line 1261 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp