OpSpringLhsBc_dU Struct Reference

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

Inheritance diagram for OpSpringLhsBc_dU:

Collaboration diagram for OpSpringLhsBc_dU:

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

Public Member Functions
	OpSpringLhsBc_dU (std::string row_field, boost::shared_ptr< SpringBcVec > bc_spring_ptr)
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)

Protected Attributes
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 OpSpringLhsBc_dU::OP = typename FormsIntegrators<FaceUserDataOperator>::Assembly<A>::OpBase

Definition at line 593 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpSpringLhsBc_dU()

OpSpringLhsBc_dU::OpSpringLhsBc_dU ( std::string  row_field, boost::shared_ptr< SpringBcVec >  bc_spring_ptr  )

inline

Definition at line 595 of file EshelbianOperators.hpp.

      : OP(row_field, row_field, OP::OPROWCOL), bcSpringPtr(bc_spring_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpSpringLhsBc_dU::iNtegrate	(	EntData &	row_data,
		EntData &	col_data
	)

Examples

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

Definition at line 2470 of file EshelbianOperators.cpp.

                                                              {
  MoFEMFunctionBegin;
 
  int row_nb_dofs = row_data.getIndices().size();
  int col_nb_dofs = col_data.getIndices().size();
  auto &locMat = OP::locMat;
  locMat.resize(row_nb_dofs, col_nb_dofs, false);
  locMat.clear();
 
  // 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()) {
 
      auto t_normal = OP::getFTensor1NormalsAtGaussPts();
      auto t_w = OP::getFTensor0IntegrationWeight();
 
      FTENSOR_INDEX(SPACE_DIM, i);
      FTENSOR_INDEX(SPACE_DIM, j);
 
      int nb_integration_pts = OP::getGaussPts().size2();
      int nb_base_functions = row_data.getN().size2();
      auto t_row_base = row_data.getFTensor0N();
 
      auto t_kd = FTensor::Kronecker_Delta<double>();
 
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
        FTensor::Tensor1<double, SPACE_DIM> t_N;
        t_N(i) = t_normal(i);
        t_N.normalize();
 
        FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_P;
        t_P(i, j) = t_N(i) * t_N(j);
        FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_Q;
        t_Q(i, j) = t_kd(i, j) - t_P(i, j);
 
        FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_d_res;
        t_d_res(i, j) = -(bc.normalStiffness * t_P(i, j) +
                          bc.tangentialStiffness * t_Q(i, j));
 
        int rr = 0;
        for (; rr != row_nb_dofs / SPACE_DIM; ++rr) {
          auto t_mat = getFTensor2FromArray<SPACE_DIM, SPACE_DIM, SPACE_DIM>(
              locMat, SPACE_DIM * rr);
          auto t_col_base = col_data.getFTensor0N(gg, 0);
          for (auto cc = 0; cc != col_nb_dofs / SPACE_DIM; ++cc) {
            t_mat(i, j) += (t_w * t_row_base * t_col_base) * t_d_res(i, j);
            ++t_mat;
            ++t_col_base;
          }
          ++t_row_base;
        }
 
        for (; rr != nb_base_functions; ++rr)
          ++t_row_base;
 
        ++t_w;
        ++t_normal;
      }
 
      locMat *= OP::getMeasure();
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcSpringPtr

boost::shared_ptr<SpringBcVec> OpSpringLhsBc_dU::bcSpringPtr

protected

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