SimpleContactProblem::OpGetContactArea Struct Reference

#include <users_modules/mortar_contact/src/SimpleContact.hpp>

Inheritance diagram for SimpleContactProblem::OpGetContactArea:

Collaboration diagram for SimpleContactProblem::OpGetContactArea:

Public Member Functions
	OpGetContactArea (const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const double cn, const bool alm_flag=false, Range post_proc_surface=Range(), double post_proc_gap_tol=0.)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
const double	cN
const bool	almFlag
VectorDouble	vecR
Range	postProcSurface
double	postProcGapTol

Detailed Description

Definition at line 3543 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetContactArea()

SimpleContactProblem::OpGetContactArea::OpGetContactArea ( const string  lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact, const double  cn, const bool  alm_flag = false, Range  post_proc_surface = Range(), double  post_proc_gap_tol = 0.  )

inline

Definition at line 3545 of file SimpleContact.hpp.

        : ContactOp(lagrange_field_name, UserDataOperator::OPCOL,
                    ContactOp::FACESLAVE),
          commonDataSimpleContact(common_data_contact), cN(cn),
          almFlag(alm_flag), postProcSurface(post_proc_surface),
          postProcGapTol(post_proc_gap_tol) {}

Member Function Documentation

doWork()

MoFEMErrorCode SimpleContactProblem::OpGetContactArea::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 4672 of file SimpleContact.cpp.

                                                                             {
  MoFEMFunctionBegin;
 
  if (data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  if (!postProcSurface.empty()) {
    const EntityHandle prism_ent = getFEEntityHandle();
    Range tri_ents;
    auto &m_field = this->getPtrFE()->mField;
    CHKERR m_field.get_moab().get_adjacencies(&prism_ent, 1, 2, false, tri_ents,
                                              moab::Interface::UNION);
    tri_ents = tri_ents.subset_by_type(MBTRI);
    if (intersect(postProcSurface, tri_ents).empty())
      MoFEMFunctionReturnHot(0);
  }
 
  const int nb_gauss_pts = data.getN().size1();
 
  int nb_base_fun_col = data.getFieldData().size();
  const double area_s =
      commonDataSimpleContact->areaSlave; // same area in master and slave
 
  vecR.resize(CommonDataSimpleContact::LAST_ELEMENT, false);
  vecR.clear();
 
  auto t_lagrange_slave =
      getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
  auto t_gap_gp = getFTensor0FromVec(*commonDataSimpleContact->gapPtr);
  auto t_w = getFTensor0IntegrationWeightSlave();
 
  for (int gg = 0; gg != nb_gauss_pts; gg++) {
    const double val_s = t_w * area_s;
    vecR[CommonDataSimpleContact::TOTAL] += val_s;
    bool gap_below_tolerance =
        postProcGapTol > std::numeric_limits<double>::epsilon() &&
        t_gap_gp < postProcGapTol;
    if (gap_below_tolerance) {
      vecR[CommonDataSimpleContact::ACTIVE] += val_s;
    } else {
      if (!almFlag &&
          SimpleContactProblem::State(cN, t_gap_gp, t_lagrange_slave)) {
        vecR[CommonDataSimpleContact::ACTIVE] += val_s;
      }
      if (almFlag &&
          SimpleContactProblem::StateALM(cN, t_gap_gp, t_lagrange_slave)) {
        vecR[CommonDataSimpleContact::ACTIVE] += val_s;
      }
    }
 
    ++t_lagrange_slave;
    ++t_gap_gp;
    ++t_w;
  } // for gauss points
 
  constexpr std::array<int, 2> indices = {
      CommonDataSimpleContact::ACTIVE,
      CommonDataSimpleContact::TOTAL,
  };
  CHKERR VecSetValues(commonDataSimpleContact->contactAreaVec, 2,
                      indices.data(), &vecR[0], ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

almFlag

const bool SimpleContactProblem::OpGetContactArea::almFlag

private

Definition at line 3561 of file SimpleContact.hpp.

cN

const double SimpleContactProblem::OpGetContactArea::cN

private

Definition at line 3560 of file SimpleContact.hpp.

commonDataSimpleContact

boost::shared_ptr<CommonDataSimpleContact> SimpleContactProblem::OpGetContactArea::commonDataSimpleContact

private

Definition at line 3559 of file SimpleContact.hpp.

postProcGapTol

double SimpleContactProblem::OpGetContactArea::postProcGapTol

private

Definition at line 3564 of file SimpleContact.hpp.

postProcSurface

Range SimpleContactProblem::OpGetContactArea::postProcSurface

private

Definition at line 3563 of file SimpleContact.hpp.

vecR

VectorDouble SimpleContactProblem::OpGetContactArea::vecR

private

Definition at line 3562 of file SimpleContact.hpp.

---

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

• users_modules/mortar_contact/src/SimpleContact.hpp
• users_modules/mortar_contact/src/impl/SimpleContact.cpp