SimpleContactProblem::OpGetGapSlave Struct Reference

Operator for the simple contact element. More...

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

Inheritance diagram for SimpleContactProblem::OpGetGapSlave:

Collaboration diagram for SimpleContactProblem::OpGetGapSlave:

Public Member Functions
	OpGetGapSlave (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	Evaluates gap function at slave face gauss points. More...

Public Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

Detailed Description

Operator for the simple contact element.

Calculates gap function at the gauss points on the slave triangle.

Definition at line 863 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetGapSlave()

SimpleContactProblem::OpGetGapSlave::OpGetGapSlave ( const string  field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact  )

inline

Definition at line 866 of file SimpleContact.hpp.

        : ContactOp(field_name, UserDataOperator::OPROW, ContactOp::FACESLAVE),
          commonDataSimpleContact(common_data_contact) {}

Member Function Documentation

doWork()

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

Evaluates gap function at slave face gauss points.

Computes gap function at slave face gauss points:

\[ g_{\textrm{n}} = - \mathbf{n}(\mathbf{x}^{(1)}) \cdot \left( \mathbf{x}^{(1)} - \mathbf{x}^{(2)} \right) \]

where \(\mathbf{n}(\mathbf{x}^{(1)})\) is the outward normal vector at the slave triangle gauss points, \(\mathbf{x}^{(1)}\) and \(\mathbf{x}^{(2)}\) are the spatial coordinates of the overlapping gauss points located at the slave and master triangles, respectively.

Definition at line 686 of file SimpleContact.cpp.

                                                                          {
  MoFEMFunctionBegin;
 
  if (data.getFieldData().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  if (type != MBVERTEX)
    MoFEMFunctionReturnHot(0);
 
  auto get_tensor_vec = [](VectorDouble &n) {
    return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
  };
 
  const int nb_gauss_pts = data.getN().size1();
 
  commonDataSimpleContact->gapPtr.get()->resize(nb_gauss_pts);
  commonDataSimpleContact->gapPtr.get()->clear();
 
  FTensor::Index<'i', 3> i;
 
  auto t_position_master_gp = getFTensor1FromMat<3>(
      *commonDataSimpleContact->positionAtGaussPtsMasterPtr);
  auto t_position_slave_gp = getFTensor1FromMat<3>(
      *commonDataSimpleContact->positionAtGaussPtsSlavePtr);
 
  auto t_gap_ptr = getFTensor0FromVec(*commonDataSimpleContact->gapPtr);
 
  auto t_normal_at_gp =
      get_tensor_vec(*(commonDataSimpleContact->normalVectorSlavePtr));
 
  for (int gg = 0; gg != nb_gauss_pts; ++gg) {
    t_gap_ptr -=
        t_normal_at_gp(i) * (t_position_slave_gp(i) - t_position_master_gp(i));
    ++t_position_slave_gp;
    ++t_position_master_gp;
    ++t_gap_ptr;
  } // for gauss points
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

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