EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS > Struct Template Reference

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

Inheritance diagram for EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >:

Collaboration diagram for EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >:

Public Types
using	OP

Public Member Functions
	OpNormalDispBcRhsImpl (std::string row_field, boost::shared_ptr< MatrixDouble > hybrid_disp_ptr, boost::shared_ptr< MatrixDouble > piola_stress_ptr, boost::shared_ptr< NormalDisplacementBcVec > &bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	hybridDispPtr
boost::shared_ptr< MatrixDouble >	piolaStressPtr
boost::shared_ptr< NormalDisplacementBcVec >	bcDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap

Detailed Description

template<AssemblyType A>
struct EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >

Definition at line 582 of file EshelbianPlasticity.hpp.

Member Typedef Documentation

OP

template<AssemblyType A>

using EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::OP

Initial value:

 
      typename FormsIntegrators<FaceUserDataOperator>::Assembly<PETSC>::OpBase

Definition at line 585 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

OpNormalDispBcRhsImpl()

template<AssemblyType A>

EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::OpNormalDispBcRhsImpl ( std::string row_field, boost::shared_ptr< MatrixDouble > hybrid_disp_ptr, boost::shared_ptr< MatrixDouble > piola_stress_ptr, boost::shared_ptr< NormalDisplacementBcVec > & bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, boost::shared_ptr< Range > ents_ptr = nullptr )

inline

Definition at line 588 of file EshelbianPlasticity.hpp.

Member Function Documentation

iNtegrate()

template<AssemblyType A>

MoFEMErrorCode OpNormalDispBcRhsImpl< A, GAUSS >::iNtegrate ( EntData & data )

protected

Definition at line 1583 of file EshelbianOperators.cpp.

                                                                       {
  MoFEMFunctionBegin;
 
  double time = OP::getFEMethod()->ts_t;
  if (EshelbianCore::dynamicRelaxation) {
    time = EshelbianCore::dynamicTime;
  }
 
  // 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()) {
 
      auto t_approx_P = getFTensor2FromMat<3, 3>(*piolaStressPtr);
      auto t_u = getFTensor1FromMat<3>(*hybridDispPtr);
      auto t_normal = OP::getFTensor1NormalsAtGaussPts();
      auto t_w = OP::getFTensor0IntegrationWeight();
 
      FTENSOR_INDEX(SPACE_DIM, i);
      FTENSOR_INDEX(SPACE_DIM, j);
 
      auto t_kd = FTensor::Kronecker_Delta<double>();
 
      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;
      }
 
      // get bc data
      double val = scale * bc.val;
 
      int nb_dofs = data.getIndices().size();
      int nb_integration_pts = OP::getGaussPts().size2();
      int nb_base_functions = data.getN().size2();
      auto t_row_base = data.getFTensor0N();
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
        FTensor::Tensor1<double, 3> t_N;
        t_N(i) = t_normal(i);
        t_N.normalize();
 
        FTensor::Tensor2<double, 3, 3> t_P;
        t_P(i, j) = t_N(i) * t_N(j);
        FTensor::Tensor2<double, 3, 3> t_Q;
        t_Q(i, j) = t_kd(i, j) - t_P(i, j);
 
        FTensor::Tensor1<double, 3> t_tracton;
        t_tracton(i) = t_approx_P(i, j) * t_N(j);
 
        FTensor::Tensor1<double, 3> t_res;
        t_res(i) = t_Q(i, j) * t_tracton(j) + t_P(i, j) * 2* t_u(j) - t_N(i) * val;
 
        auto t_nf = getFTensor1FromPtr<3>(&*OP::locF.begin());
        int bb = 0;
        for (; bb != nb_dofs / SPACE_DIM; ++bb) {
          t_nf(i) += (t_w * t_row_base) * t_res(i);
          ++t_nf;
          ++t_row_base;
        }
        for (; bb != nb_base_functions; ++bb)
          ++t_row_base;
 
        ++t_w;
        ++t_normal;
        ++t_u;
        ++t_approx_P;
      }
 
      OP::locF *= OP::getMeasure();
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcDispPtr

template<AssemblyType A>

boost::shared_ptr<NormalDisplacementBcVec> EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::bcDispPtr

protected

Definition at line 603 of file EshelbianPlasticity.hpp.

hybridDispPtr

template<AssemblyType A>

boost::shared_ptr<MatrixDouble> EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::hybridDispPtr

protected

Definition at line 601 of file EshelbianPlasticity.hpp.

piolaStressPtr

template<AssemblyType A>

boost::shared_ptr<MatrixDouble> EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::piolaStressPtr

protected

Definition at line 602 of file EshelbianPlasticity.hpp.

scalingMethodsMap

template<AssemblyType A>

std::map<std::string, boost::shared_ptr<ScalingMethod> > EshelbianPlasticity::OpNormalDispBcRhsImpl< A, GAUSS >::scalingMethodsMap

protected

Definition at line 604 of file EshelbianPlasticity.hpp.

---

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

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