OpBrokenPressureBc Struct Reference

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

Inheritance diagram for OpBrokenPressureBc:

Collaboration diagram for OpBrokenPressureBc:

Public Member Functions
	OpBrokenPressureBc (const std::string row_field, boost::shared_ptr< PressureBcVec > bc_data, boost::shared_ptr< double > piola_scale_ptr, boost::shared_ptr< MatrixDouble > hybrid_grad_disp, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv)
MoFEMErrorCode	iNtegrate (EntData &data)

Protected Attributes
boost::shared_ptr< PressureBcVec >	bcData
boost::shared_ptr< double >	piolaScalePtr
boost::shared_ptr< MatrixDouble >	hybridGradDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap

Detailed Description

Examples

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

Definition at line 499 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpBrokenPressureBc()

OpBrokenPressureBc::OpBrokenPressureBc ( const std::string  row_field, boost::shared_ptr< PressureBcVec >  bc_data, boost::shared_ptr< double >  piola_scale_ptr, boost::shared_ptr< MatrixDouble >  hybrid_grad_disp, std::map< std::string, boost::shared_ptr< ScalingMethod > >  smv  )

inline

Definition at line 501 of file EshelbianOperators.hpp.

      : FormsIntegrators<FaceUserDataOperator>::Assembly<PETSC>::OpBase(
            row_field, row_field, FaceUserDataOperator::OPROW),
        bcData(bc_data), piolaScalePtr(), hybridGradDispPtr(hybrid_grad_disp),
        scalingMethodsMap(smv) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpBrokenPressureBc::iNtegrate

(

EntData &

data

)

Examples

mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1969 of file EshelbianOperators.cpp.

                                                          {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(3, i);
 
  int nb_dofs = data.getFieldData().size();
  int nb_integration_pts = getGaussPts().size2();
  int nb_base_functions = data.getN().size2();
 
  double time = getFEMethod()->ts_t;
  if (EshelbianCore::dynamicRelaxation) {
    time = EshelbianCore::dynamicTime;
  }
 
#ifndef NDEBUG
  if (this->locF.size() != nb_dofs)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Size of locF %ld != nb_dofs %d", this->locF.size(), nb_dofs);
#endif // NDEBUG
 
  auto integrate_rhs = [&](auto &bc, auto calc_tau, double time_scale) {
    MoFEMFunctionBegin;
 
    auto val = bc.val;
    auto t_row_base = data.getFTensor0N();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    auto t_tangent1 = getFTensor1Tangent1AtGaussPts();
    auto t_tangent2 = getFTensor1Tangent2AtGaussPts();
 
    auto t_grad_gamma_u = getFTensor2FromMat<3, 2>(*hybridGradDispPtr);
 
    double scale = (piolaScalePtr) ? 1. / (*piolaScalePtr) : 1.0;
 
    for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
      FTENSOR_INDEX(SPACE_DIM, i);
      FTENSOR_INDEX(SPACE_DIM, j);
      FTENSOR_INDEX(SPACE_DIM, k);
      FTensor::Number<0> N0;
      FTensor::Number<1> N1;
 
      FTensor::Tensor1<double, 3> t_normal;
      if (EshelbianCore::stretchSelector == LINEAR &&
          EshelbianCore::gradApproximator < MODERATE_ROT) {
        
        t_normal(i) = (FTensor::levi_civita<double>(i, j, k) * t_tangent1(j)) *
                      t_tangent2(k);
      } else {
        t_normal(i) = (FTensor::levi_civita<double>(i, j, k) *
                       (t_tangent1(j) + t_grad_gamma_u(j, N0))) *
                       (t_tangent2(k) + t_grad_gamma_u(k, N1));
      }
      auto tau = calc_tau(t_coords(0), t_coords(1), t_coords(2));
      auto t_val = FTensor::Tensor1<double, 3>();
      t_val(i) = (time_scale * t_w * tau * scale * val) * t_normal(i);
 
      auto t_f = getFTensor1FromPtr<3>(&*this->locF.begin());
      int rr = 0;
      for (; rr != nb_dofs / SPACE_DIM; ++rr) {
        t_f(i) += t_row_base * t_val(i);
        ++t_row_base;
        ++t_f;
      }
 
      for (; rr != nb_base_functions; ++rr)
        ++t_row_base;
      ++t_w;
      ++t_coords;
      ++t_tangent1;
      ++t_tangent2;
      ++t_grad_gamma_u;
    }
    this->locF /= 2.;
 
    MoFEMFunctionReturn(0);
  };
 
  // get entity of face
  EntityHandle fe_ent = getFEEntityHandle();
  for (auto &bc : *(bcData)) {
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
 
      double time_scale = 1;
      if (scalingMethodsMap.find(bc.blockName) != scalingMethodsMap.end()) {
        time_scale *= scalingMethodsMap.at(bc.blockName)->getScale(time);
      }
 
      int nb_dofs = data.getFieldData().size();
      if (nb_dofs) {
        CHKERR integrate_rhs(
            bc, [](double, double, double) { return 1; }, time_scale);
      }
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcData

boost::shared_ptr<PressureBcVec> OpBrokenPressureBc::bcData

protected

Definition at line 514 of file EshelbianOperators.hpp.

hybridGradDispPtr

boost::shared_ptr<MatrixDouble> OpBrokenPressureBc::hybridGradDispPtr

protected

Definition at line 516 of file EshelbianOperators.hpp.

piolaScalePtr

boost::shared_ptr<double> OpBrokenPressureBc::piolaScalePtr

protected

Definition at line 515 of file EshelbianOperators.hpp.

scalingMethodsMap

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

protected

Definition at line 517 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