SimpleContactProblem::OpGetNormalSlave Struct Reference

Computes, for reference configuration, normal to slave face that is common to all gauss points. More...

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

Inheritance diagram for SimpleContactProblem::OpGetNormalSlave:

Collaboration diagram for SimpleContactProblem::OpGetNormalSlave:

Public Member Functions
	OpGetNormalSlave (const string field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	Evaluates unit normal vector to the slave surface vector based on reference base coordinates. More...

Public Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact

Detailed Description

Computes, for reference configuration, normal to slave face that is common to all gauss points.

Definition at line 664 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetNormalSlave()

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

inline

Definition at line 667 of file SimpleContact.hpp.

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

Member Function Documentation

doWork()

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

Evaluates unit normal vector to the slave surface vector based on reference base coordinates.

Computes normal vector based on reference base coordinates based on mesh (moab vertices) coordinates:

\[ {\mathbf N}^{(1)}({\boldsymbol{\chi}}(\xi, \eta)) = \frac{\partial\mathbf{X}(\xi, \eta)}{\partial\xi}\times\frac{\partial \mathbf{X}(\xi, \eta)} {\partial\eta} \]

where \({\boldsymbol{\chi}}^{(1)}(\xi, \eta)\) is the vector of reference coordinates at the gauss point on slave surface with parent coordinates \(\xi\) and \(\eta\) evaluated according to

\[ {\boldsymbol{\chi}}(\xi, \eta) = \sum\limits^{3}_{i = 1} N_i(\xi, \eta){{\boldsymbol{\chi}}}_i \]

where \( N_i \) is the shape function corresponding to the \( i-{\rm{th}}\) degree of freedom in the reference configuration \({{\boldsymbol{\chi}}}^{(1)}_i\) corresponding to the 3 nodes of the triangular slave face.

Definition at line 384 of file SimpleContact.cpp.

                                                                             {
  MoFEMFunctionBegin;
 
  if (data.getFieldData().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  if (type != MBVERTEX)
    MoFEMFunctionReturnHot(0);
 
  FTensor::Index<'i', 3> i;
 
  auto get_tensor_vec = [](VectorDouble &n) {
    return FTensor::Tensor1<double *, 3>(&n(0), &n(1), &n(2));
  };
 
  const double *normal_slave_ptr = &getNormalSlave()[0];
 
  commonDataSimpleContact->normalVectorSlavePtr.get()->resize(3);
  commonDataSimpleContact->normalVectorSlavePtr.get()->clear();
 
  auto t_normal =
      get_tensor_vec(*(commonDataSimpleContact->normalVectorSlavePtr));
 
  for (int ii = 0; ii != 3; ++ii)
    t_normal(ii) = normal_slave_ptr[ii];
 
  const double normal_length = sqrt(t_normal(i) * t_normal(i));
  t_normal(i) = t_normal(i) / normal_length;
 
  commonDataSimpleContact->areaSlave = 0.5 * normal_length;
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataSimpleContact

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

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