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

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

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

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

Public Types
using	OP

Public Member Functions
	OpBrokenPressureBcLhsImpl_dU (std::string row_field, boost::shared_ptr< PressureBcVec > bc_data, boost::shared_ptr< MatrixDouble > hybrid_grad_disp, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)

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

Detailed Description

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

Definition at line 531 of file EshelbianPlasticity.hpp.

Member Typedef Documentation

OP

template<AssemblyType A>

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

Initial value:

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

Definition at line 534 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

OpBrokenPressureBcLhsImpl_dU()

template<AssemblyType A>

EshelbianPlasticity::OpBrokenPressureBcLhsImpl_dU< A, GAUSS >::OpBrokenPressureBcLhsImpl_dU ( std::string row_field, boost::shared_ptr< PressureBcVec > bc_data, boost::shared_ptr< MatrixDouble > hybrid_grad_disp, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv )

inline

Definition at line 537 of file EshelbianPlasticity.hpp.

                      {
  ExternalStrain(std::string name, std::vector<double> attr, Range ents);
  std::string blockName;

Member Function Documentation

iNtegrate()

template<AssemblyType A>

MoFEMErrorCode OpBrokenPressureBcLhsImpl_dU< A, GAUSS >::iNtegrate ( EntData & row_data, EntData & col_data )

protected

Examples

EshelbianOperators.cpp.

Definition at line 2091 of file EshelbianOperators.cpp.

                                                                     {
  MoFEMFunctionBegin;
 
  if (EshelbianCore::stretchSelector == LINEAR &&
      EshelbianCore::gradApproximator < MODERATE_ROT) {
    MoFEMFunctionReturnHot(0);
  }
 
  double time = OP::getFEMethod()->ts_t;
  if (EshelbianCore::dynamicRelaxation) {
    time = EshelbianCore::dynamicTime;
  }
 
  int nb_base_functions = row_data.getN().size2();
  int row_nb_dofs = row_data.getIndices().size();
  int col_nb_dofs = col_data.getIndices().size();
  int nb_integration_pts = OP::getGaussPts().size2();
  auto &locMat = OP::locMat;
  locMat.resize(row_nb_dofs, col_nb_dofs, false);
  locMat.clear();
 
auto integrate_lhs = [&](auto &bc, auto calc_tau, double time_scale) {
    MoFEMFunctionBegin;
 
    auto val = bc.val;
    auto t_row_base = row_data.getFTensor0N();
    auto t_w = OP::getFTensor0IntegrationWeight();
    auto t_coords = OP::getFTensor1CoordsAtGaussPts();
    auto t_tangent1 = OP::getFTensor1Tangent1AtGaussPts();
    auto t_tangent2 = OP::getFTensor1Tangent2AtGaussPts();
 
    auto t_grad_gamma_u = getFTensor2FromMat<3, 2>(*hybridGradDispPtr);
    auto t_kd = FTensor::Kronecker_Delta<double>();
 
    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_INDEX(SPACE_DIM, l);
 
      FTensor::Number<0> N0;
      FTensor::Number<1> N1;
 
      auto tau = calc_tau(t_coords(0), t_coords(1), t_coords(2));
      auto t_val = time_scale * t_w * tau * val;
 
      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_diff_col_base = col_data.getFTensor1DiffN<2>(gg, 0);
        for (auto cc = 0; cc != col_nb_dofs / SPACE_DIM; ++cc) {
          FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_normal_du;
          t_normal_du(i, l) = (FTensor::levi_civita<double>(i, j, k) *
                           (t_tangent2(k) + t_grad_gamma_u(k, N1))) *
                              t_kd(j, l) * t_diff_col_base(N0)
 
                          +
 
                          (FTensor::levi_civita<double>(i, j, k) *
                           (t_tangent1(j) + t_grad_gamma_u(j, N0))) *
                              t_kd(k, l) * t_diff_col_base(N1);
 
          t_mat(i, j) += (t_w * t_row_base) * t_val * t_normal_du(i, j);
          ++t_mat;
          ++t_diff_col_base;
        }
        ++t_row_base;
      }
 
      for (; rr != nb_base_functions; ++rr)
        ++t_row_base;
      ++t_w;
      ++t_coords;
      ++t_tangent1;
      ++t_tangent2;
      ++t_grad_gamma_u;
    }
 
    OP::locMat /= 2.;
 
    MoFEMFunctionReturn(0);
  };
 
  // get entity of face
  EntityHandle fe_ent = OP::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 = row_data.getFieldData().size();
      if (nb_dofs) {
        CHKERR integrate_lhs(
            bc, [](double, double, double) { return 1; }, time_scale);
      }
    }
  }
 
  MoFEMFunctionReturn(0);
 
}

Member Data Documentation

bcData

template<AssemblyType A>

boost::shared_ptr<PressureBcVec> EshelbianPlasticity::OpBrokenPressureBcLhsImpl_dU< A, GAUSS >::bcData

protected

Definition at line 546 of file EshelbianPlasticity.hpp.

hybridGradDispPtr

template<AssemblyType A>

boost::shared_ptr<MatrixDouble> EshelbianPlasticity::OpBrokenPressureBcLhsImpl_dU< A, GAUSS >::hybridGradDispPtr

protected

Definition at line 547 of file EshelbianPlasticity.hpp.

scalingMethodsMap

template<AssemblyType A>

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

protected

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