OpSpringRhsBc Struct Reference

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

Inheritance diagram for OpSpringRhsBc:

Collaboration diagram for OpSpringRhsBc:

Public Types
using	OP = typename FormsIntegrators< FaceUserDataOperator >::Assembly< PETSC >::OpBase

Public Member Functions
	OpSpringRhsBc (std::string row_field, boost::shared_ptr< MatrixDouble > hybrid_disp_ptr, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< SpringBcVec > bc_spring_ptr)
MoFEMErrorCode	iNtegrate (EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	hybridDispPtr
boost::shared_ptr< std::vector< BrokenBaseSideData > >	brokenBaseSideDataPtr
boost::shared_ptr< SpringBcVec >	bcSpringPtr

Detailed Description

Examples

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

Definition at line 590 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpSpringRhsBc::OP = typename FormsIntegrators<FaceUserDataOperator>::Assembly<PETSC>::OpBase

Definition at line 593 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpSpringRhsBc()

OpSpringRhsBc::OpSpringRhsBc ( std::string  row_field, boost::shared_ptr< MatrixDouble >  hybrid_disp_ptr, boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_side_data_ptr, boost::shared_ptr< SpringBcVec >  bc_spring_ptr  )

inline

Definition at line 596 of file EshelbianOperators.hpp.

      : OP(row_field, row_field, OP::OPROW), hybridDispPtr(hybrid_disp_ptr),
        brokenBaseSideDataPtr(broken_side_data_ptr), bcSpringPtr(bc_spring_ptr) {
  }

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpSpringRhsBc::iNtegrate

(

EntData &

data

)

Examples

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

Definition at line 2338 of file EshelbianOperators.cpp.

                                                     {
  MoFEMFunctionBegin;
 
  // get entity of face
  EntityHandle fe_ent = OP::getFEEntityHandle();
  // iterate over all boundary data
  for (auto &bc : (*bcSpringPtr)) {
    // check if finite element entity is part of boundary condition
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
      
      for (auto &bd : (*brokenBaseSideDataPtr)) {
 
        auto t_approx_P = getFTensor2FromMat<3, 3, -1, DL>(bd.getFlux());
        auto t_u = getFTensor1FromMat<3, -1, DL>(*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>();
 
        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_traction;
          t_traction(i) = t_approx_P(i, j) * t_N(j);
 
          FTensor::Tensor1<double, 3> t_res;
          t_res(i) = 0.5 *(t_traction(i)) - bc.normalStiffness * t_P(i, j) * t_u(j) -
                     bc.tangentialStiffness * t_Q(i, j) * t_u(j);
 
          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 * OP::getMeasure()) * 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;
        }
      }
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcSpringPtr

boost::shared_ptr<SpringBcVec> OpSpringRhsBc::bcSpringPtr

protected

Definition at line 609 of file EshelbianOperators.hpp.

brokenBaseSideDataPtr

boost::shared_ptr<std::vector<BrokenBaseSideData> > OpSpringRhsBc::brokenBaseSideDataPtr

protected

Definition at line 608 of file EshelbianOperators.hpp.

hybridDispPtr

boost::shared_ptr<MatrixDouble> OpSpringRhsBc::hybridDispPtr

protected

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